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"Beam me over," Alice: A cricket's quantum journey 

This thesis addresses two known quantities in quantum information science: 
(1) entanglement cost, and (2) Holevo capacity. These quantities will be 
crucial values when teleportation becomes common in daily life, perhaps 
centuries from now. 

Assume that Alice desires to send a singing Japanese cricket to her friend 
Bob in America, and that Alice and Bob already share a quantum entan- 
glement. First, Alice sends Bob a mass of information bits resulting from 
the interaction between the cricket she holds in her hand and half of the 
entanglement. Subsequently, Bob receives the information bits and ma- 
nipulates the other half of the entanglement, transforming them back into 
the original cricket. Examining this situation from an instrumental engi- 
neering viewpoint, quantifying the amount of the quantum entanglement 
and the number of information bits is crucial for this transmission. If both 
values are enough, Alice could even send herself to Bob's place instead of 
the tiny cricket. 

The topics of this thesis therefore are: (1 ) the mathematical properties of 
the entanglement cost, such as whether it is an additive measure similar to 
normal length or weight; and (2) how to calculate the Holevo capacity, an 
ultimately achievable limit of the information conveyance capacity of an 
information channel, such as of a single photon passing through an optical 
fiber or space. These two distinct quantities are magically tied together by 
several "additive or not" hypotheses, which await mathematical proof. 
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ABSTRACT 



Quantum entanglement is essential in quantum information science. Quantum computations require quantum 
entanglement, and that is utilized in various quantum communication protocols. This dissertation focuses on the 
theory of quantum entanglement, for bipiirtitc mixed states in terms of the asymptotic behavior, ranging from 
quEuitification of the entanglement to the Holevo capacity of quantum channels. 

Quantifying quantum entanglement for a given quantum bipartite state p can be done by two ways: one 
is by counting the number of the Bell pairs required to produce p, and the other is by counting the number 
of the Bell pairs that can be extracted from p. This dissertation focuses primarily on the values given by the 
former way, which is called the entanglement cost denoted by Ec- Calculating Ec is doubly difficult; Ec is 
represented in the limitation form Ec{p) = lim„_,(x) -EF(p®")/n, while Ep being called the entanglement of 
formation appearing in the previous formula is represented in an optimization form on a high-dimensional and 
seemingly multimodgil function. To overcome these difficulties, the author sought whether Ep is additive, ngimely, 
Ef{pi ® P2) = Ef{pi) + Ef{p2) holds or not for bipartite quantum states. If this equality holds for every case, 
Ef = Ec holds, which makes the problem much feasible. 

In this dissertation, the value Ef and its additivity for antisymmetric states on a 3®3-level system are 
discussed in order to calculate Ec- This is to find the first example to be calculated the value of Ec for 
nontrivial certain states since none of nontrivial examples were known to be calculated Ec when this resesirch 
was performed. Furthermore, the inequality of the strong superadditivity of Ef and the calculation of the Holevo 
capacity are discussed because the inequsJity holding, the additivity holding for the Holevo capacity, and the 
additivity holding for Ef are closely related. Trials for finding the counterexamples are numerically performed for 
the inequalities in the case of the nontrivially lowest-dimensional space. This space is composed of twin bipartite 
two-level systems, which is consequently I6-dimensional. 

Finally, calculating the Holevo capacity on qubit chginnels is performed, seemingly for the first time. The 
original formula of the Holevo capacity is an optimization form of a multimodal function, tiius iiard to calculate 
straightforwardly. The presented algorithm solves that optimization with a fixed arbitrary small error, maximizing 
a high-dimensional concave function corresponding to the specified error. The relation between the computation 
cost and the error is analyzed. 
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Prologue 
— Introduction and Background — 




Chapter 1 



Quantum Information Science 



Quantum information science, the consolidation of the rules and the intelligence of this physical world, 
is the emerging essence of science. We need tools to understand and investigate this developing science. 

Prior to the chapter of the introduction, we present here what the quantum information science 
is in this chapter. Concepts used in this dissertation are explained, compactly, neatly, and hopefully 
coherently. 

1.1 What is a quantum state? 

What is a quantum state? It specifies a specific physical state at a given time in the framework 
of quantum physics. Quantum physics attempts to explain most everything in the world, especially the 
micro-world. Generally, a quantum state is the "superposition" of multiple specified states, as often 
depicted by Schrodinger's cat, that is the superposition of a living cat and a dead cat. One needs 
to understand such peculiar phenomena because information devices are becoming smaller and smaller 
toward nanometer world. Predictions say that fifteen or twenty years from now, each bit of information 
will be contained in as few as one atom as technology continues to grow exponentially, as it has for more 
than forty years. 

In this dissertation, we deal solely with static quantum states; we do not consider physical con- 
tinuous time transition on states that might involve the Schrodinger equation, that is, H{t)\^{t)) = 



\/—lh-g^\'ip{t)) . Still, in the static framework, there are many interesting phenomena such as quantum 

entanglement. It is an anomalous correlation on two sites or more, which is not depicted in the framework 
of classical views of physics. Einstein et al. opposed quantum physics because quantum entanglement 



ment is, however, being experimentally confirmed, and many protocols utilizing it are proposed and 
being tested experimentally. 

Here, again, we ask "What is a quantum state?" The definitions, mainly from mathematical view- 
points, follows. 




seemed to raise 



the issue of super luminal communication jEmstem-Podoisky- Rosen 35|_ Quantum entangle- 
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1.1.1 Pure state and its tensor product 
Pure states 

A pure state \^) is a vector of which the length is one, dweUing on a specified vector space, or a 
Hilbert space, of a complex number field. One can represent it by a column vector of which the elements 
are complex numbers as, 

/ -1 \ 



X2 



with the conditions xi,X2, ■ ■ ■ ,Xd e C and \xi\'^ + |a;2p + ... + \xd\'^ = 1. Note that this is just a 
representation, and a vector treated in quantum information science is an element of a vector space, 
more precisely a Hilbert space over complex number field. The space is spanned by some bases, such as 
{1^); {IT)i U)}i or simply |2), . . . , as exemplified in the next. 

The polarization of a single photon: The polarization of a photon flying in z-direction is conceptually 
I I 

y 



polarizing 
plate 




Figure 1.1: This figure depicts how to prepare photons with the polarization with cos 9 |^) + sin 6 ||). 

represented in a two-dimensional space with bases |<-!-) and ||), that represent z-directional and y- 
directional polarization, respectively. If the polarization angle is 9 on the xy-plane, the state of the 
photon is represented as \9) — cos 9 |^)+sin0 ||) that is the superposition oi |<-!-) and ||) with coefficients 
cos 9 and sin 9. This vector \9) is represented as (sf,jg). See Fig ll.ll Note that the coefficients can be 
any complex number as long as their squared summation is one. cos 9 |^) + e^^^^ sin 9 ||) is such an 
example, which can be considered to have the phase shift of 7 in the factor of |, and it is called to have 
circular or elliptical polarization. 

As a physical aspect, the dimensionality of two is enough for the vector space of the polarization 
of a photon. In principle, the most possible measurement of the polarization of each single photon is 
equivalent to check whether the photon passes through a prepared polarizing plate, and the passed 
photon loses the information of the polarization. Thus all the information the observer can get about 
the polarization of the photon is as few as two possibilities, which is related to the dimensionality of two 
here. 

The magnetic moment of a single silver atom: Each silver atom has its momentum state on the 
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y spin up 



S 



t) 



silver beam 




furnace 



N 



spin down 



i) 



Figure 1.2: Stern-Gerlach experiment 



vector space spanned by jt) and The state is written as a\1) + f3\i) that is a superposition of jf) 
and II). Physically, the momentum can be measured by making the atom flying through nonuniform 
magnetic field, and the atoms curve in two possible directions. The measurement is done by observing 
which of the two directions the atom has curved. If the magnetic field is set up to detect whether the 
atom has the moment of |t) or of |J,), then the atom turns to become |t) or ||) by curving into the 
corresponding directions with the possibility \a\'^ or respectively. This experiment is called Stern- 
Gerlach experiment , which is a typical experiment to show the physical phenomenon of quantization of 
the magnetic moment of atoms. 

A physical system is called a d-level system when the system is a space for d-dimensional vectors. 

The physical system to be observed the polarization of a single photon or the magnetic momentum 
of a single silver atom is a typical 2-level system. 

Ket and bra 

As a mathematical convenience, \ip) is called a ket vector. {ip\ is the Hermitian transpose, or the 
conjugate transpose of lip), which is called a bra vector. The former is considered a column vector, and 
the latter is considered a row vector on matrix arithmetic. Thus, one can consider to be a square 

matrix. 

Tensor product of pure states 

You might like to consider the state of a system which contains two particles and more. Here, we 
give the concepts of tensor product for two pure states. The tensor product H^H' of two spaces H and 
H' is the d x d'-dimensional vector space when H and H' are d- and d'- dimensional, respectively. The 
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tensor product of two vectors of H and H' is defined as follows: 



For example, ^ 



X2 



2/2 



\ xd J \ yd' J 



( xiyi \ 
xiyi 



xiVd' 
X2yi 

X2yd' 

\ XdVd' ) 



(1.1) 



/0.4 




By convention, \<^\) |02) is sometimes abbreviated 



as or <^2). One can consider that tensor product of states of two particles is considered to be 

the state of the whole two particles. 



1.1.2 Mixed states and tensor product 

A mixed state is represented by a semi-positive Hermitian matrix^ of which the trace, or the 
summation of the diagonal elements, is 1. It is regarded as a stochastic mixture of multiple pure states 
{\tlJi)}i, as p = J2iPi\'^i){'^i\ with {pi}i a probability distribution. Note that stochastic mixture is a 
different notion from superposition of pure states. The semi-positive Hermitian matrix of a mixed state 
is called the density matrix. This dissertation employs the conventional phrase "a state p on H" or "a 
state p on C^" if p is a d x d semi-positive matrix of which the trace is one. Notation B{'H) is employed 
to explicitly specify the set of states on a Hilbert space H. 

Notation 1.1 (mixed states) 

B{'H) is the set of any mixed states p on a Hilbert space H. When H is d-dimensional, B{H) is equal to 
the set of dx d semi-positive matrix of which the trace is one. 



Note that a pure state \<p) is regarded to be equal to a mixed state as a stochastic mixture of 

a single state with 100% weight. Mixed states other than pure states may be called non-pure states. 
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The tensor product of mixed states is an expansion to the tensor product of pure states, as follows: 
(ail ■■■ aid\ ( hii ■■■ bid' \ 



( aiibii ■ ■ ■ aiibid' 



aiibd'i ■ ■ ■ cLiibd'd' 



adibii ■ ■ ■ adibid' 



\ adib, 



d'l 



^dlbd'd' 



\ adi ... add J \ bd'i ■ ■ ■ bd'd' J 



aidbii 



addbi 



aidbd'i ■ ■ ■ aidbd', 



addb 



Id' 



addbd'i ■ ■ ■ addbd'd' j 



(1.2) 



1.1.3 Arithmetic on states 

This subsection deals with some arithmetic on quantum states, such as, tracing out, von Neumann 
entropy and quantum divergence. 

Notation 1.2 (trace out, tracing out) 

For a state p = J2iPi pf ® pf on Ti-A ^Ti-B, tracing out operations are defined as follows. Tracing out 
the space of Hb from p is defined as 

TrBp = ^p,(Trpf).pf, (1.3) 

i 

resultingly a state on TCa- Similarly, tracing out the space ofTCA is, 

Tr^p = ^p,(Trpf).pf on Hb- (1.4) 

i 

Ttb is an operation tracing out the space of B, leaving the space of A. Conversely, Tia is an 
operation tracing out the space of A, leaving the space of B. Concrete examples of tracing out appear 
at ppinOl 

Notation 1.3 (logp, y^) 

For a semi-positive Hermite matrix p, log p and ^fp are deRned as follows: 



log p^U diag [log ti, . . . ,\ogtd]U^ 

for p = f/diag [ti, . . . ,td] U\ 

and 

VP = C/ diag [^i,...,^U^ 

for p^U diag [ii, . . . , td] U\ 

where diag [ti, . . . , td] is the diagonal matrix ■ . , and U is a unitary matrix. 



(1.5) 



(1.6) 
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To keep consistency and to avoid confusion the following conventions are employed here: 

logO is treated as — oo, additionally, OlogO is treated as 0. (1.7) 

Notation 1.4 (von Neumann entropy, j Meumann 32 ]^ 

S{p)^-Trp\ogp (1.8) 

One can show that S{p) is equal to the Shannon entropy of the whole eigenvalues of p, that is, 
S{p) — log Ai) with {Xi} being the eigenvalues of p. Note that for pure states \(j)){(p\, the von 

Neumann entropy is zero. 

Notation 1.5 (quantum divergence) 

The quantum divergence for two mixed states from the same space are defined as follows: 

Hip\\a)=TTiplogp-plogcT). (1.9) 

Quantum divergence is also called quantum relative entropy. Quantum divergence has following 
properties, as it is similar to a measure of distance of two states, but it lacks the property of symmetry. 

Proposition 1.6 (basic properties of i/( || )) 



H{p\\a)>0, (1.10) 

Hip\\a)^0 ^ p = a, (1.11) 

H{;\\-2)^H{-2\\-i). (1.12) 

1.1.4 Mapping on quantum states 

This subsection deals with mapping on quantum states. In this dissertation, a "map" or "mapping" 
means a linear map. One can associate a mapping as a time transition of a quantum state, formalized 
as a quantum channel later. 

Notation 1.7 (tensor product of mapping) 

For two mappings Ai and A2, the tensor product Ai (X) A2 is defined based on a restriction that Ai (Xi 
A2(pi P2) — Ai(pi) (g) A2(/02) for any input pi of Ai and any input p2 of A2. 

Notation 1.8 (complete positivity, CP-ness) 

A complete positive map A is a map such that for any (finite) dimensional identity map /, A (g) / maps a 
semi-positive Hermitian matrix into a semi-positive Hermitian matrix. This complete positivity is also 
called CP-ness. 



Notation 1.9 (trace preserving, TP-ness) 

A trace preserving map A is a map satisfying Tr Ap = 
also called TP-ness. 



Tr p for any p. This trace preserving property is 



Notation 1.10 (CPTP-ness) 

Combining the two terms above, a map of CP and TP is a CPTP map. 
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Notation 1.11 (Bures distance and the fidelity) 

The Gdelity between given two states is deGned as 



F(p,a) = Tr^VP^VP- (1-13) 

The Bures distance is dehned as 



B{p,c7) ^ 2^1 -F{p, a). (1.14) 

Proposition 1.12 

The fidelity has the properties such as: 

F{-i,-2) < 1 (1.15) 

F(-i,-i) = 1 (1.16) 

F{;,-2) = F{-2,-i) (1.17) 

F{\^){^\,\(l)){(l)\) = for the same dimensional (1.18) 

The fidelity is considered to be how truly a quantum state is transmitted. 100% fidelity, or the fidelity 
being one is the perfect transmission, which means the transmission is done without losing information 
the quantum states held. Note that in l|1.18(l . the left hand side is determined even if the dimension of 
IV') and 10) is different while the right hand side is not defined. 

The Bures distance (I1.14f) is a derivation from the fidelity to take on the properties of a distance. 
The nearer upward to 1 or 100% the fidelity between two states become, the smaller downward to the 
distance between the two states become. 



1.2 A quantum channel and its capacity 

There are various frameworks to deal with time transition of a quantum states as follows: 



(1) Schrodinger equation — 




(2) Quantum circuit — 




H 




and, 




(3) Quantum channel — 


- 


\ v 

channel • 



(1) and (2) treat only pure states, and (3) treats mixed states generally. To investigate quantum 
information science, treating mixed states is necessary. 

Definition 1.13 (quantum channel) 

A quantum channel A is a CPTP map (see Notation M.li^) . Thus, (Ap) ® a is a quantum state for any 
quantum state a. (See Fis \l.t!l ) 

Because of the linearity of mappings of quantum channel, p' — A(p) is abbreviated as p' — Ap. 
Definition 1.14 (Tensor product of channels) 

For channels A and A', the tensor product of these channels A (g) A' is a linear map which maps p® p' to 
(Ap) (g) (A'p') where p is any input of A and p' is any input of A'. 

The concept of tensor product is natural to consider multiple channels in the physical world. 



8 



-.0 .■ channel O 



O ^ O 

before after 

Figure 1.3: A quantum channel must be a CPTP map in order to keep the whole world a quantum state 
that is represented as a semi-positive Hermitian matrix. 

Theorem 1.15 (Stinespring's dilation) 

For any quantum channel A : p i— s- p', tiiere exists some auxiliary space K , some auxiliary state a, and a 
unitary matrix U such that 

p' = TrKU (p®(j)U\ (1.19) 

Now we have formalized what a quantum channel is. Then, how is the capacity of the channel 
considered, i.e., how much information can be carried from the sender to the receiver in a remote place 
through a quantum channel? One of the formalizations is the Holevo capacity. The definition follows. 

Definition 1.16 (Holevo capacity, j Ho'evo 73 | |Holevo 98 ||Sch umacher- Westmoreland 97 |) 

The Holevo capacity C of a given quantum channel A is as follows: 

n n 

C(A)= max SiY p., Ap^) - Y S{Ap,), (1.20) 

i=l 1=1 

wiiere rt G > 0, X^ILi^'* ~ ^ every pi s is an input of the channel A. 

The optimized n in H1.20|l is regarded as the number of kinds of the input states of the channel A, 
when the communication capacity attains the Holevo capacity of A in an asymptotic sense. 

Theorem 1.17 (i"hya-^etz-Watanabe 9f|^ 

C(A) = minmaxi/(Ap|| Act) (1.21) 

a p 

where <j and p are the inputs of the channel A. 

Physically, the Holevo capacity is a classical information capacity of a given quantum channel at 
which input particles are not allowed to be entangled with each other, and the output particles are 
measured collectively. To fully utilize quantum aspects of a quantum channel, one might consider the 
capacity at which the input particles are freely entangled with each other to send a message. This 
capacity is called the full capacity. 

Definition 1.18 (Full capacity) 

C'(A) = lim C(A®")/n (1.22) 
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1.3 Entanglement and its quantification 

This section defines the quantum entanglement in the case of a bipartite system. 

Wlicn considering bipartite quantum entanglement, we often think of two figures, named Ahce and 
Bob. Their spaces, at which quantum states are considered, are represented as Hilbert spaces Ha and 
Hb, respectively. 

Definition 1.19 (entanglement) 

For a pure state \(p) £ Ha ^ Hb, |0) is said to be entangled 

if 10) cannot be represented in a form ® \(j)B) with £ Ha and |(/)_b) £ Hb- 
For a mixed state p on Ha ® Hb , P is said to be entangled 
if p cannot be represented as J2iPi\'f'i)i't'i\ 

with uncntangled pure states and a probability distribution {pi}. 

An unentangled state is called separable. To explicitly refer to two system which are entangled, notation 
Ha '■ Hb is used in this dissertation. When considering the entanglement, our interest reside at two sites 
(Alice and Bob) rather than their space Ha and Hb- In such a case we just denote Alice: Bob or A:B. 

Definition 1.20 (LOCC, Cf. jBennet-DiVicenzo-Smolin-Wootters 9 6| ^ 

LOCC operations (standing for Local Operation plus Classical Communication) for a bipartite system 
are dehned to be comprised only of: 

(1 ) local CPTP operations on each side, and 

(2) communications between the two sites by only classical means. 

Classical communication means the following: One side physically measures the state of its own system, 
and transmit the output of the measurement to the other side. Based on this information, the receiver 
may perform physical operations of its own system. Physical measurements are formalized to be the 
concept of POVM, which we do not take up in this dissertation. 




Figure 1.4: LOCC operations: Only Local Operations and Classical Communications are allowed. Joint 
operations or sending quantum information is prohibited if it is beyond the allowed operations. 

There is a well-established entanglement measure for bipartite pure states. 

Notation 1.21 (reduced von Neumann entropy) 

The reduced von Neumann entropy for a bipartite state is defined as 

E''--^{m = s{TTB\m) 

for a pure state £ Ha ® Hb- (1-23) 
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One can easily conErm that S{Ttb \(f>M) = 5'(Tr^ \(f>){(j>\) thus E '^'■^ = S'(Tr^ 

If a pure state is entangled, the reduced von Neumann entropy is larger than 0. 
Example 1.22 (Bell state, Bell pair) 

are called the Bell states. These states |<i>^), 1$^), l^*^) and l^*") are often considered to be implicitly 
shared by Alice and Bob throughout this dissertation. Their reduced von Neumann entropy is 1. 

Bell states are often considered to be units of quantum entanglement. 

For bipartite mixed states, there are various candidates to measure quantum entanglement as follows. 

Definition 1.23 (entanglement cost, entanglement distillation i Bennet-DiVicenzo-Smoiin- Wootters 96 | ^ 

The entanglement cost is defined as 

E^--\p) = inf { e I ^(e, 5)^ (m, n, L) : |e - ^| < 5, B (i(|vI/-)(*-|^™), < e } , (1.25) 

and the entanglement distillation is defined as 

E^-^{p) - sup { e I ^(e, 5)^ (m, n, L) : |e - ^ | < 5, B (|*-)(vi/- L(p«")) < e } , (1.26) 
where e > 0, e, (5 > 0, m, n e N and L{-) is a LOCC operation. B{-, •) is the Bures distance defined in 

n~m . 

The entanglement distillation Ed is the asymptotic quantity of Bell states distilled from p with 
LOCC operations. Thus, one can consider Eo as a measure to quantify the entanglement as a resource 
to be used for quantum teleportation or quantum super dense coding. 

The entanglement cost Ec is the asymptotic quantity of Bell states necessary to produce p with 
LOCC operations. 




Distillation 

Figure 1.5: Difference between Ec and Ed 



Definition 1.24 (entanglement of formation i Bennet-DiVicenzo-Smoiin- Wootters 96 | ) 

Tiie entanglement of formation is defined as 

n 

EP-^{P)^ min ^Ki?'^^^(|0,)(0,|), (1.27) 

n; \(pl),...,\4)r^}■,Pl,■■■,Pr^ ^ 

with p, > 0,J2"=lPi = '^^T,^=lP^\<P^){<P^\ = P- 
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We may use the denotations E{-), Ec{-), Ei:,{-), Ep{-) which omit the superscription A : B when the 
two sites are clear to consider entanglement. 



Proposition 1.25 ( iHayden-Horodecki-TerhaToTl ^ 

The following holds: 

Ec{p) = lim 



(1.28) 



This equality H1.28|l is significant in that the conceptually defined Ec represented as H1.25|l has 
become mathematically defined by substituting H1.27|l into H1.28|l . The calculation is, however, not 
simple, which leads to the main subject of this dissertation because (ll.28|) being substituted with (|1.27() 
is in a limitation form over infinitely many optimization forms {i?F(p'*")/«}$^Li- 



1.4 Example: a two- level system 

This section presents properties of the Bloch sphere representing states on a quhit, or a 2-level 
system.^ 

1.4.1 Bloch sphere and the Stokes parameterization 

Assume 



p{x,y,z) = 



(1.29) 



1 / 1 + z X — V — ly 
X + V— Ty 1 — z 

is a 2-dimensional mixed state. One can easily confirm that the semi-positivity is equivalent to x'^ + 
+ z"^ < 1. The representation of a 2-dimensional mixed state under this condition is called the Stokes 
parameterization. The sphere {(x, y, z) \ x'^ + y^ + z'^ < 1}, in which each point is associated with the 
mixed state p{x,y,z), is called the Bloch sphere. The surface of the Bloch sphere corresponds to pure 
states. 




Figure 1.6: The polarizations of a photon on the Bloch sphere 
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1.4.2 A qubit channel 

A qubit channel, which maps a qubit space to a qubit space, can be regarded as an afhne trans- 
formation on the Bloch sphere. Thus a qubit channel is characterized by the output ellipsoid inside the 
Bloch sphere. See Fig. 11.71 
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Figure 1.7: A qubit channel's output ellipsoid 



1.4.3 The von Neumann entropy on the Bloch sphere 

In order to consider the Holevo capacity of a qubit channel, here we consider the von Neumann 
entropy on the Bloch sphere. The eigenvalues of p{x, y, z) fsee ll.29|l are ''"^"^^ and , 

so the von Neumann entropy is the binary entropy of them. The von Neumann entropy on the Bloch 
sphere is spherically-symmetric and concave function. 

1.5 Relation between Ep and the Holevo capacity 

Due to Stinespring's dilation jstinespring 55| ^^^^y CPTP map A : ^(Hin) ^ ^(Hout) can be represented 
as the composition of an isometric embedding-^ of B{Hin) into a bipartite system B{Hout) S('^aux) 
followed by the operation tracing out S(7Yaux) leaving the output space i3(7iout), represented as, 

A: (1.30) 

p ^ UpW ^ Trn^^^UpU^ 

by choosing Tiaux of which the dimension large enough and the isometric embedding U . 
Notation 1.26 {K-x, image vectors of the isometric embedding) 

Denote K-a ~ UHin, a subspace ofHout "Waux, the image subspace ofU where U appeared in H1.3Up . 

We can say that the channel A is equivalent to a tracing out operation from JC\ to the output space 
'Hout, with an isometric embedding U. 
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Figure 1.8: Left: The von Neumann entropy on the xy-section of the Bloch sphere. Note that the 
domain of this function is originally three-dimensional, and it is reduced to two-dimensional xy-plane in 
this figure. Right: The binary entropy. 



Theorem 1.27 (see [Matsumoto-Shimono- winter 041 ^ 

C{A) = sup{5(Tr„_p) - : p on JCa}- (1.31) 

where C(-) is the Holevo capacity defined in ^1.21]) . 

Notation 1.28 {0\{c B{ICx)), preimage matrices of the average output.) 

For a given channel, denote as follows: 

Oa = argmax{5(Trw_/5) - £;r"^^"^(/9) : p on ICa}. (1.32) 

pee(/CA) 

Theorem 1.29 

If C(A (g) A') = C(A) + C(A') holds, then following holds: 

Vpe C'a,Vp' e Oa' £;r'^^"^(/5® p') = ^^r'^^'^^'lp) + ^r'^^"''(p')- (1-33) 

□ 

It is interesting whether C(A'*") = nC{A) holds for a given/every channel A, because one can 
conclude whether Epip'^"') = nEp{p) holds by the theorem above. As partial results, following channels 
are known to satisfy this additivity property: 

• Unital qubit channels" , cf. j^ing 02c||gTng02a| 

• Entanglement -breaking channels^ , cf. j ^hor 02 | 

• Arbitrary depolarizing channels King^ , cf. jKmg 02b| 
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1.6 Additivity problems 



This section gives issues related to the additivity of Ep and to the Holevo capacity. 
Definition 1.30 (strong superadditivity of Ep) 

For a state p on Hai ® 'Ha2 ®'Hbi ®'Hb2, the inequality of strong superadditivity of Ep is defined as, 

EP--^{P) > EP-^{Tt2 p) + EP-^{Tt^ p). (1.34) 

Here Tr2 means tracing out the space ofTiAi ^Ti-Bi 8,nd Ti'i means tracing out the space ofTi.A2 ^ 'Hb2- 
See Fig. 11. .91 Note that whether this inequahty holds for every case has not yet solved. 




Figure 1.9: Illustration of the strong superadditivity. Entanglement is measured (left) Al + A2 : Bl + B2, 
(right) Al : Bl and A2 : B2 separately. The inequality (11.3411 means the quantity of the left is larger than 
the sum of the quantities of the right. 

Theorem 1 31 ( jlVIatsumoto-Shimono-Winter 041 IShor 041 ^ 

Tiie following four propositions are equivalent: 

• Tiie additivity in Ep : ^{p, a) : Ep{p®a)^ Ep{p) + Ep{a) 

• Tiie additivity in the Holevo capacity : ^(Ai, A2) : C(Ai ® A2) = C(Ai) + C(A2) 

• Tiie strong superadditivity of Ep for any mixed states 

• The strong superadditivity of Ep for any pure states 

1.7 Miscellaneous 

Here, we fix some conventions, which might otherwise cause confusion. 

1.7.1 Base of logarithm 

Throughout this dissertation, the base of logarithm is fixed to two (2) unless the base is specifically 
indicated. Therefore, log 2 = 1. Note that this rule is applied to the logarithm of matrices like logp as 
well, which effects von Neumann entropy and the quantum divergence and so on. 

1.7.2 Convex and concave 

Convexity of a function: a function / is convex when it bends downward, as 

/(Ax + (1-A)y) < A/(x) + (l-A)/(y) for < A < 1. (1.35) 
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Concavity of a function: a function / is concave when it bends upward, as 

/(Ax + (1 - X)y) > Xfix) + (1 - X)f{y) for < A < 1. 



(1.36) 




convex concave 
Figure 1.10: A convex function (Left) and a concave function (Right) 

1.7.3 Tensor products 

In this chapter, or in this dissertation, the tensor product is defined for 

• pure states, 

• mixed states, 

• mapping or quantum channels 

1.7.4 The terms level and dimension 

In quantum informatics, the terms level and dimension may not be different in nature as one can 
say a state comes from 2-diniensional space or 2-level system. In this dissertation, however, they may 
be distinguished in some context, as the level refers to the quantum space and the dimension refers to 
technical derivation to perform the optimizing calculation for the Holevo capacity. 

Endnotes of Chapter [T] 

1 Semi-positive Hermitian matrix: A matrix p is Hermitian as long as it satisfies p = . In this 
case, all its eigenvalues are real numbers. It is semi-positive when all its eigenvalues are equal to or larger 
than zero. 

2 "Qubit" can be considered as a unit of quantum information as the counterpart of the concept 
of "bit" of the current information theory. Qubit is named by Benjamin Schumacher. 

3 A linear operator or map U is called an isometric embedding if [/^[/ is an identity operator. The 
name isometric embedding emerged from the following two reasons: 

• Any vector |0) does not change its length after the mapping [/, as || \(t>) \\ = \\ U\(f>) \\. This is 
because |i |0) f - = m^U\(j)) = (C/|</.))t([/|0)) = || um \\' and || |0) || > 0, || U\^) \\ > 0. 
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• Any vectors do not change their distance after the mapping U as || |0) — \tjj) \\ = \\ U\(f)) — U\i(}) \\. 

Note that isometric embedding is not always a unitary transformation, even though any unitary trans- 
formation is always an isometric embedding. 

4 A unital channel is a quantum channel A : B(Win) — > ^(Wout) which maps the maximally mixed 
state of Win to the maximally mixed state of Wout- The maximally mixed state is defined as ^ Ym=\ I*X^I 
for BiC^). 

5 An entanglement-breaking channel A : yB(7iin) B{'Hout) is an entanglement breaking channel if 
it satisfies the following condition. Any output of A (g) 1 : B{Hin) (Xi B{Hs.nx) B{Hout) B{Hs.nx} is not 
entangled between Tiout and TYaux for any auxiliary space Tiaux- Here 1 is the identity map of fi(7Yaux)- 

6 A channel A : B{TL) B{Ti) is an depolarizing channel if it is a unital channel and there exist 
s € R such that A((t) — A(p) = s {a — p), namely, A is a proportionally shrinking mapping toward the 
maximally mixed state. 
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ket 
bra 

B{T-L) the set of mixed states of a given Hilbert space 
transpose of a given matrix 

conjugate transpose or Hermitian transpose of a given matrix 
-y, -z) Stokes parameterization in the Bloch sphere 



/9, a : mixed quantum states 
B{p, a) Bures distance between two states 
F{p, a) fidelity between two states 

H[p II cr) quantum divergence, quantum relative entropy 
S{p) von Neumann entropy of a state 



•, □ : a mixed quantum state share by the two sites A and B 
E^'^{-) reduced von Neumann entropy 
E^-^{-) entanglement cost 
^D^i') entanglement distillation 
Ep'^{-) entanglement of formation 
^N^i') logarithm of negativity 
Er^{') relative entropy of entanglement 

partial transpose of a given matrix 



A : a quantum channel 

C(A) Holevo capacity of a quantum channel 
C{A) full capacity of a quantum channel 

/Ca image vectors induced by an isometric embedding ?7 of a channel A 
Oa preimage matrices of the average output of A in B{JCk) 



Table 1.1: Denotations used in this dissertation. 
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Chapter 2 



Introduction to the dissertation 



The mam body of this dissertation is the author's researches in quantum information science with a main 
interest in quantum entanglement, and methodologies of conducting research as well. 

Before entering the chapters of the research performed, we are going to depict (i) the background 
of this dissertation, (ii) how the research was conducted, and (iii) this dissertation's organization. 



2.1 Background 

2.1.1 Quantum information science 

Quantum information science, the consolidation of information science and quantum physics, both of 
which date back to the first half of the twentieth century, has attracted the interest of many researchers 
over the last ten years as the invention of efficient algorithms for number factorization^ [$hor 94 ] ^ 
database search^ j crover 96 | ^ j'I^qj-q jg ^ characteristic phenomena in quantum physics that is called quan- 
tum entanglement, which plays an essential role in quantum computation jjozsa-Linden 02 |^ Various quan- 
tum protocols were devised using quantum entanglement, such as quantum teleportation'^ IsBgjP'vy ^ 
and quantum superdense coding* Bennett- wiesner 92 |^ 

The author is concerned with the entanglement cost iBennet-DiVicenzo-Smoiin-Wootters 96 | ^ ^^^^ ,^1^^ 
Holevo capacity i Hoievo 73 | ^ j^-^^ entanglement cost is one way to quantify the entanglement of a bipartite 
quantum state, and the Holevo capacity is a classical communication capacity of a given quantum 
channel. These two different measures are related by j Matsumoto-Shimono-Winter o D through Stinespring's 
dilation theorem as presented in section [TTSl 



2.1.2 Quantifying entanglement 

How to quantify the entanglement of quantum states? Quantum entanglement raise up useful effects 
on quantum computation/communication mysteriously, thus quantification of entanglement is significant 
and natural methodology to clarify the mystery. 

To measure the usefulness of entanglement, there is a simple way is to compare two entangled 
pairs. There may be a situation that a pair particles shared by Alice and Bob is more entangled than 
another pair particles also shared by the two figures because the former pair is more useful to be used in 
quantum communication, say, sending a message in quantum superdense coding. Then how to compare? 
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An established way is to see the possibility of transformation from a pair to the other through LOCC 
operation, as LOCC cannot create additional entanglement. If a pair #1 can be transformed into another 
pair #2 by LOCC, then one can say that #1 is more useful than #2. 

jVidai 99 1 g[jQ-^^g good theorem for this question for pure states. Two pairs of pure states can be 
compared by the concept of majorizatiort' This majorization was the first significant step to measure 
quantum entanglement. 

To quantify quantum entanglement literally in quantity, an asymptotic comparison is conceived. 
The idea is, for a state p, to investigate the possibility to transform from/into by LOCC 

with large N and M. |v]>~)'**^ is considered to have M ebits. The idea is, if p^^^ can be transformed 
into Ivt/^)*^*^^ then p is more than Ni/Mi ebits, and if p^^^ can be transformed from then p 

is less than N2/M2 ebits. The former idea is formalized in entanglement distillation (Ejj) and the latter 
is formalized in entanglement cost {Eq) (see H1.25|l and (|1.26(l ). 




Distillation 



Figure 2.1: The depiction of the entanglement cost and the entanglement of distillation. The transfor- 
mations are LOCC. p is compared with copies of a Bell state, say, l^"") — . 

For pure states iP°pescu-Rohriich 97 | ividai o ol iNieisen ooi|Donaid-HorQdecki-Rudoiph 02| ^^le entanglement cost 
and the entanglement distillation coincide each other. For mixed states, they are, however, different. 
That is, irreversibility occurs in asymptotic transformation through LOCC for mixed bipartite state. 
Thus quantifying entanglement is not simple for mixed states. Another idea is to put a restrict such as a 
measure should be additive. There has been devised many measures and the entanglement of formation 
(Ep) is also a good candidate to measure entanglement. It is originally defined as a formula H1.27|l . and 
the relation (|1.28|) connects it to the entanglement cost. Even though the entanglement distillation seems 
to break the additivity |Shor-Smoiin-Terhai m \ ^ ^j^g entanglement cost and the entanglement of formation 
seem not break the additivity, or it is hard to find such a state p that breaks the additivity. 

2.1.3 Channel capacity beyond classical physics 

For a classical channel, that is a stochastic mapping from the input alphabets to the output alpha- 
bets, there is a sole capacity formalized by Shannon, which gives the supremum quantity of information 
that is transmitted asymptotically per channel use EUaiiOll. To determine the capacity for classical 
case, there are efficient algorithms ^M^sMM^M, 

For a quantum channel, there could be many capacities for a channel. An interesting quantity is the 
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Holevo capacity, which was upper-bounded by the Holevo quantity j Hoievo 73 |^ ^^^^ proved to be equal to 

the Holevo quantity [Holevo 981 ISchumacher- Westmoreland 971 

There were proposals to calculate the Holevo capacity for general quantum channels josawa-Nagaoka oi | ^hor 
jOsawa-Nagaoka oi| expauds the method of the Arimoto-Blahut algorithm. 1 ^^™ made a comment utiliz- 
ing the column generation method. This dissertation presents seemingly the first algorithm to calculate 

the Holevo capacity of qubit channels jHayashi-lmai-Matsumoto-Huskai-Shimono 04| _ 

The Holevo capacity is a measure which does not allow entanglement among input particles. There 
is a possibility that allowing entanglement among the input particles might increase the capacity, even 
though the Shannon capacity keeps the additivity. This problem is formalized as whether C(A'^") — 
n C(A) for n = 2, 3, . . . , or C(Ai (g) A2) = C(Ai) + C(A2) holds for arbitrary channels. This additivity 
was related to the additivity of Ep i Matsumoto-^himono-winter g H^ Verifying the additivity of the Holevo 
capacity is an interesting problem, thus numerically calculating the Holevo capacity for high-dimensional 
case is expected. 



2.2 Research flow of this dissertation 

This dissertation is concerned with the problems of quantification of quantum entanglement and 
the computation of the Holevo capacity. As is written above, quantifying the entanglement of bipartite 
pure states is almost settled, in that the reduced von Neumann entropy is the sole measure with the 
appropriate conditions. The appropriate properties here are the reducing property against any non- 
quantum physical operation (formalized in LOCC), additivity property, and so on. 

2.2.1 Quantifying Entanglement 

Quantum states are generally mixed states, as each of them presented in the stochastic mixture of 
some pure states. Various ways of quantifying the entanglement of a given mixed state are proposed 
iB(^nn(^t-DiVicenzo-Smoiin-Wootters 9^ 1 ^ ^^^^i as the entanglement cost (Ec), the entanglement of formation 
(Ep), the entanglement distillation {Ed), and the relative entropy of entanglement [E^)^ . For a given 
p, Ec{p) is the measure of how many the Bell states are needed to construct p in the asymptotic sense. ^ 
Ep is similarly defined, but rather technically. Ep) of p is defined as how many the Bell states can be 
distilled from p. 

Naive questions might arise: "Which is the best [in quantifying the entanglement of bipartite mixed 
states]?", "What is the appropriateness in measures?", "Is there only one proper measure?" 
A relation 

Ec{p) = lim (2.1) 



n — *oo 



,s given by jHayden-Horodecki-ierhaJoij ^ ^^^^^ ^j^g following expectation were brought about: 

• Ec could be calculated analytically and/or numerically. 

• The additivity of Ep, that is 

Ep{p'^^)=' nEp{p), (2.2) 
might holds thus Eq = Ep holds and Ec is calculated rather easily. 
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If the additivity does not hold, calculating Ec is a formidable problem. Calculating 



Ef{p) 



mm 

n;pi,...,p„;|</)i>,...,|0„) 



i=l 



P^E{\c^,){<i>,\) 



n e N, 

Pl, . . . ,Pn > 0, E"=lP» = 1. 

|0i), . . . , |(/)„) are pure states, 



(2.3) 



where i? is the reduced von Neumann entropy, is quite an high-dimensional optimization problem, even 
if the dimension of p is quite low such as two or three^ . Still more, Eq is a limitation form of this 
problem. 

The first trial of calculating Ec for nontrivial case was done for three-level antisymmetric states 
[Vidai-Diir-girac Qi j _ -y^ith the expectation of calculating the first case of Ec for nontrivial case, the author 
gave lower bounds of antisymmetric states Ishimono 2] ^ .^^j^j^-jj appear in Chapter 3. More general form 
of lE2J) is 

EF{p(?)cr) =■ EFip)+EF{a). (2.4) 

The author seem to gave the first specific example j shimono o al ^ presented in Chapter 4. Conclusively, the 
first examples to calculate Ec for nontrivial case was done by i vidai-Diir-Cirac 02 |^ ^^j. -j-j^g mixture of Bell 
states. Ec for three-dimensional antisymmetric states was calculated by E^^^ and the result was 
further expanded to more general antisymmetric case j Mats^moto-Yura q4 | _ 

As of the days when Ec and Ed were proposed i Bennet-Diyicenzo-gmoiin-Wootters 9^ 1 ^ whether they can 
be different or not was unknown^ . Such states were found by j Horodecki ^ j^j, g ^ 3-level mixed states. 
Attempts to find other examples was done by iV'dai-Dur-Cirac 02 | ^ -pj^g g^^j-j^Q], ^^j-^jg dissertation also gave 
examples utilizing the following facts: 

• Ef can be calculated through the Holevo capacity iMatsumoto-ghimono-Winter q I^ 

• Ed is bounded by another easily-calculable quantity called the logarithmic negativity [Y'dai-Wemer 92 1 
-Eat, i.e. Ed < En 

This result is presented in Chapter 5. 



2.2.2 Strong superadditivity 

There is another quantification called "the relative entropy of entanglement" jvedrai-Pienio-Kippin-Knight 97| 
denoted En. It was believed to be additive for a while, however, negated later by finding counterexamples 

jVollbrecht- Werner oTl 

It seems that determining only one measure is not appropriate because Ec{p) > Ed{p) holds for 
some state p (examples is presented in Chapter 5), even though both Ec and Ed are good candidates to 
measure quantum entanglement. The additivity is the property E{p ® p') — E{p) + E{p') or E{p^") — 
n E{p). It is of natural and expected requirement for any kind of "measuring" . As quantum entanglement 
is the valuable peculiarity from an ordinal classical physical viewpoint, one might like to quantify it like 
counting money. From a theoretical viewpoint of convenience, if the additivity holds for Ep, then 
Ep = Ec, which leads to a difficult formula of Ec in calculation to a more feasible one calculating Ep. 

Based on the concept written above, the value of Ep and Ec for some specific states are tried 
calculated. For example, the antisymmetric states of bipartite three-level systems are taken up. (These 
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states can be regarded as two particles of fermions in one system, thus consideration of this entanglement 
might cause practical applications.) In order to try to calculate Ec{p) where p is an antisymmetric state, 
Ef{p'^'^) is calculated. (This work is completed by 1 ^"'''^ as Ec{p) = 1, which means Ef{p'^^) for 
Vn e N is calculated. Furthermore, it is extended to the space of ® {£d^^m-2 ^^^^SSmEM^ as 
the value is log2 d — \ ) 

With the theorem of Stinespring dilation, the article of j Matsumoto-Shimono-Winter 4] related the Holevo 
capacity to Ep. Utilizing this fact, the author calculated Ep for a 2® A - level system, which shows gaps 
between Ep and Ep,. To check the additivity of Ep in the lowest nontrivial case, the inequality of the 
strong superadditivity for a twin 2 ® 2-level system^" . Feasibility in calculation of the von Neumann 
entropy for a 4-level system and calculation of Ep for 2 2 j Hiii-Wootters 9 7|| Wootters 98 ] jg utilized, even 
though calculating Ep seems formidable even for a 2^ (g) 2^-level system. However, the inequality of the 
strong superadditivity is easy to check rather than the additivity check, thus the attempts to check this 
is performed in Chapter 

2.2.3 The Holevo capacity 

As the Stinespring correspondence introduces, the action of any quantum channels (CPTP transfor- 
mation) is equivalent to the result of the following procedurel ^'^*""'""*"-^^™"'^"-^'^*'^'- g^ l (ggg Yig. 1^ : 

- making a tensor product from the input state and a certain auxiliary state, 

- transforming it by a certain unitary operator, and then 

- reducing into the space of the output. 

Ep of the intermediate state is equal to j Matsitmotp-ghirnQno-Wint^^r q D. 

the von Neumann entropy of the average output minus the Holevo capacity. (2.5) 

Therefore, if we can calculate the Holevo capacity, we can calculate the Ep for various states. In 
this dissertation, given examples are calculated utilizing a mathematical programming package NUOPT 
[nvqpt i q£ Mathematical Systems Inc., which is capable of solving large optimization problems, and the 
convergence in calculation is analyzed. 

outer environment £ 



input p intermediate ct output hp-Jro 

Figure 2.2: This figure shows the Stinespring correspondense. Any quantum channel can be represented 
as releasing some space £ after a unitary operation over the input state and an additional space. 

The Holevo capacity is a classical information capacity of a given memoryless quantum channel. 
Back to the history, it was formalized and upper-bounded by the Holevo quantity 1^°^*^™ . The Holcvo 
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capacity was proved to be equal to the Holevo quantity later j Hoievo 98 |Eih umacher-Westmoreiand 97 |^ ^j^j^j^ 
is defined as, 

n n 

C{A)= max 5( Vp, Aft) - Vp, 5(ApO (2.6) 

for a quantum channel A : B{TL±-a) B{Ti.ont), where S{-) is the von Neumann entropy. 
In 1998, another representation was provided as johya-Petz-Watanabe 97]. 

C(A) = min max H(Ap\\Aa). (2.7) 

where -ff(-||-) is the quantum divergence. 

Difficulty in calculation 

Both of (|2.()|l and (|2.7|) contain the optimization operators such as min and max, and the functions 
to be optimized seem challenging to optimize because: 
- The objective function H2.6|l is convex w.r.t. {pi)i. 

Therefore, it is hard to guarantee that a found local maximum is the global one. 

There are, however, feasible ways to calculate the Holevo capacity numerically.^^ The objective 
function of H2.6|) is concave w.r.t. eachpj because the first term is concave as S is concave and the second 
term is linear w.r.t. each pi. Therefore, once every pi is fixed, it is easy to find the global maximum by 
gradual descending. The author employs the methodology of covering with a fine lattice the whole space 
to specify {pi)i beforehand for qubit channels (see Fig l2.3|l . The finer the lattice becomes, the closer the 
approximative solution becomes to the actual value. 



1 -iT^^A/-*^ 

as y c , » t , . V» 
• • » • ^ . ' » » * 


^^^^^^^^^^^^ 
Jf * * * t « e ? n a -J ^ ' 4 ' ' 
'"it t ^ ^ t ^ * " i * * J.*!* 
I- t t iL J J 




\ \t ^ I I * * ' 

"D \ ^^^-"^ 
\ ^^.^^'^^ 05 

° \ ..-^^''"^'^^ 
V-"''''''^ -0.5 


'»jsA ' t * " 3 3 J i */J 

. - - 5 

.^-^ 
-0,5 \^^.^^--^-o.Ei 

i40 





Figure 2.3: The Bloch spheres covered with lattices. Each equator and meridian of them are divided 
into 20, 40 and 100 sections respectively. Specifying {pi)i as the vertices of a lattice beforehand, 1)2.61) is 
feasible to solve because it becomes a problem of maximizing a concave function with a convex search 
space w.r.t. {pi)i- The finer the lattice becomes, the closer the solution with the lattice becomes to the 
actual Holevo capacity. 
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Results 



The results in terms of the Holevo capacity of this dissertation are as foUows: 

- The Holevo capacity of qubit channels became possible to calculate 31 years after the formula 
was proposed. 

- S{Y:tlP^ M) - E7=lP^ S{Ap^) is maximized by NUOPT iHnnPll. 

- Qubit channels which require four inputs to achieve their Holevo capacity are found. 

- The convergence of the optimization calculation is analyzed for 

- the algorithm maximizing S{J2'i=iPi ^Pi) - J27=iP^ S{Api). 

2.2.4 Additivity issues 

Here in this subsection, the additivity issues related with quantum entanglement is described. 

iMatsumoto-ghimono-Winter Q4| j-evealcd the relation between Ep (the entanglement of formation) and 
the Holevo capacity which shows the additivity dependence between them. This relations became to be 
called the MSW correspondence, and using this relation, 1^^°'' found additivity equivalence among 
various problems, including the strong superadditivity of Ep. 

As aforementined, additivity of measures is important; additivity means that the quantity measured 
for two objects equals the sum of the quantities for the two measured separately, and this additivity 
property depends on the measure how to quantify. Back to the topic, whether Ep is additive or not 
is important because the additivity of Ep leads the additivity of Ec and no one knows whether Ep is 
additive or not. Additionary, other measures of entanglement such as Ep) (entanglement distillation) 
and Efi (the relative entropy of entanglement) were shown to be non-additive, and researchers may like 
to find out the measure for quantum entanglement that satisfies the additivity. Parallel to this, they like 
to find out the capacity of quantum channel that satisfies the additivity, and the Holevo capacity is the 
strong candidate, even though entanglement between input particles may increase the communication 
efficienty more than the sum of the communication efficiency using non-entangled particles, but nobody 
knows whether this capacity holds or not. 

The importance of the equivalence properties iM ^tsumoto-Shimono-Winter 04 |[ shor 04 | jg follows. 

// the additivity of Ep holds: the strong superadditivity of Ep holds, and the additivity of 
the Holevo capacity holds. 

// the additivty of Ep does not hold, which means there some two quantum states p and a 
satisfying Ep{p (g) cr) < Ep{p) + Ep{a): the strong superadditivyt of Ep breaks for some 
quantum states, and the additivity of the Holevo capacity breaks (for some tow channels Ai 
and A2, C(Ai (g) A2) > C(Ai) + C(A2) with the aide of entanglement between input particles 
of the two channels). 

From the viewpoint of this additivity, this dissertation can be seen as follows: Chapter 3 and 4 tried 
to check whether Ep is additive for the limited, highly symmetric states. Chapter 6 and 7 tried to find 
the witness against the additivty, through trying to find the counterexample against the superadditivity 
of Ep and trying to find the breakage of the Holevo capacity of a certain unordinary channel, respectively. 
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2.3 Organization of this dissertation 

The organization of this dissertation is as foUows: 

Part I consists of Chapter 1 and this Chapter 2, mentioning the dissertation's introductory 
matters. In Chapter 1, we review the basic notations and the definitions in quantum infor- 
mation theory to be used in subsequent chapters. This Chapter 2 is the introduction for the 
following main body of this dissertation. 

Part II, consisting of Chapters 3, 4, 5 and 6, deals with the topics related to entanglement 
quantification and its additivity problems. Chapters 3 and 4 deal with the entanglement cost 
of antisymmetric states. Chapter 5 deals with the gap between Ec and Ed. Chapter 6 deals 
with the strong superadditivity. 

Part III, consisting of Chapter 7, deals with the topics of calculating the Holevo capacity. The 
Holevo capacity is a significant quantity defined more than thirty years ago, the numerical 
calculation is, however, difficult. We present the algorithms to calculate it and application of 
it. 

Part IV, the final part, consist of Chapter 8, summarize this dissertation. 

Plenty parts of this dissertation are derivation from journal articles and conference presentations, 
such as : 

Chapter 3 from Chapter 4 from j shimono 03 |^ Chapter 5 from j Matsumoto-Shimono-vyinter o D^ 

Chapter 7 from i^^y^^^^~^^^^~^^^^^"^Q^'-'~^^^^^~^^^™°^Q Q^l 



Endnotes of Chapter [H 

1 Shor's algorithm : it solves the factorization in O(n^) with n being the number of digits, even though 
the fastest classical algorithm does it in approximately 0{exp{n^^^)). Note that the former is polynomial and 
the latter is superpolynomial. 

2 Grover's algorithm: it finds the target in a database in size n with the time proportional to 0(n^^^). 
This means the time cost grows only 10 times if the database size grows as much as 100 times. 

3 Quantum teleportation: it teleports quantum states to a remote place with classical communication and 
quantum entanglement between the two sites beforehand. 

4 

Quantum superdense coding: it transmits 2 bits while sending 1 qubit with prepared quantum entanglement. 
5 Majorization: two pairs of pure states l^pi) and |i/)2) are compared as follows: 

- 1. Let {Xi} and {ni} be each of the eigenvalues in descending order of Tr^ l^iXV-"!! and TrA \i'2){i'2\, 
respectively. 

- 1. Compare the sequences {Ai, Ai + A2, Ai + A2 + A3, . . .} and {/^i, ^1 + M2, A*i + + MSi ■ • ■}• 
They are equivalent 

- when can be transformed to 1-02) by LOCC and 

- when all of Ai < /^i, Ai + A2 < /ii + /i2, Ai + A2 + A3 < /^i + M2 + /is, . • . hold. 
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This method to compare is called majorization. 

6 The relative entropy of entanglement Er: Er[p) is defined as min H{p\\a). 

fT: separable 

7 The asymptotic sense here means that p is not single to be compared with the Bell states. Rather, p®" 
for large n is compared with m copies of the Bell states, and ^ is considered in the limit of n, m — > oo. 

8 For a bipartite d-dimensional system, n in the formula 12.31 could reach d? , and each of {4>i}i^i is a 
d^-dimensional vector on the bipartite system on C, represented in 2ci^ ~ 1 parameters in R, as well as {p^jf^i 
being considered — l)-dimensional parameters in R. One can remove one parameter of a global phase of a 
quantum system. As a whole, the number of parameters is (2d^ — 1) x d'^ + [d? — 1) — 1 = 2d'' — 2. Then, how 
to optimize such a high-dimensional problems!? It seems very difficult to calculate Bp even for d = 2. (See the 
table below.) Fortunately, HiU-Woottors 97 Woottors 98 ^^^^ analytical formula for d = 2 to calculate Ep. 



d 


(1) 2 


3 


4 


5 6 


7 


2d'^-2 


(0) 30 


160 


510 


1248 2590 


4800 . . . 



Table 2.1: The dimensional explosion in computing Ep 



9 The difference between Ec{p) and Ed{p) means some LOCC operations are irreversible (even) in an 
asymptotic sense. Namely, bulk of the Bell state in quantity Ec{p) may be transformed into bulk of p in 

unit quantity by some LOCC operations, but any LOCC operation cannot retrieve the original quantity of the 
Bell state (refer to the definition of Ec and Eo-) 

10 

- Firstly, to consider quantum entanglement for a bipartite system, 2 ® 2 is the lowest-dimensional case; 
for an n ® m-level system, n < 2 or m < 2 is trivial. 

- Secondly, to check the additivity, twofold (i.e. ®^) of 2 ® 2 is the minimum. 

Therefore, to check the additivity, the nontrivial lowest case to calculate the quantity of entanglement is (2(8)2)®^ 
dimensional, which becomes 16 dimensional. 
11 

The calculation was done in 2004 using an existing calculator; quantum computers are not yet available! 
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Figure 2.4: The contributions of this dissertation. 
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Part II 

Entanglement 
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Problems concerning quantification of quantum entanglement are integrated in this part. The first 
two chapters deal with the problem of the entanglement of formation of antisymmetric states in a 3 (E) 3- 
level system. Next, irreversibility of entanglement is shown through presenting the difference in Ep 
and Ed for a mixed 2 (g) 4-level system. Lastly, the strong superadditivity for a (2 (g) 2)®^-level system 
is numerically verified, which gave circumstantial evidences for the additivity of the entanglement of 
formation. 

Chapters in Part II 

3 Vp on : Ec{p) > 0.585 , or Epip'^") > 0.585?^. 

- Attempts for calculating the Ec- 

4 Wp onHl : Epip) ^ 1,Ef{p'^^) ^2. 

- Checking the additivity in Ep. 

5 3p on <C?®C^: Ed{p) < Ec(,p). 

-~ The gap shows the irreversibility in LOCC. 

6 3p on (C2 ® C2)^' : Ep{p) ^ Ep{1:ti p) + Ep{1:t2 p)? 

- Attempts for seeking any counterexample against the strong superadditivity 
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Chapter 3 



Lower Bounds for Entanglement Cost of 
Antisymmetric States 

Calculating Ec is quite difficult. A breakthrough has been desired. We paid attention to antisymmetric 
states. Lower bounds were found even though the exact values were unknown. 

This chapter gives a lower bound of the entanglement cost for antisymmetric states of a bipartite 
d-level system to be log2 ^jrj . This chapter is a derivation from I shimono 2] ^ 

3.1 First step toward calculating Ec 

The entanglement cost can be determined by asymptotic behavior of the entanglement formation 
|Hayden-Horodecki-Terhaioi|^ regarded to bc very difficult to calculate. 

This chapter gives a lower bound of the entanglement cost of antisymmetric states for a bipartite 
d-level system. We prove that all of the eigenvalues of each reduced matrix of any pure state affiliating to 
7i®" for general d is not greater than (^x'")"' '^here H- is the set of antisymmetric states of a bipartite 
d-level systems, defined in the following section. This is proved by investigating a CP map A defined in 
(HI. 

3.2 Problem setup 

The author makes the following assumptions: Ha and TCb are d-dimensional Hilbert spaces with 
the basis D = {\i)}i=i...d- Hab — Ha He- For 1 < i < j < d, 

= e Uab ■ 

D' — j))}i<i<^j<d- The antisymmetric space = span_D' C Hab ■ 
Notation 3.1 (matrices) 

For a positive integer m, 9Jl,„ is a set of m x m-dimensional matrices with each entry a complex num- 
ber C. For a set X, [ayljjgjf is a matrix of which any (i, j) -component specified aij, and dJl{X) = 
{['^ij]i,jex\{o-ij} C C} is a collection of matrices with each row and column labeled with elements of X. 
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Notation 3.2 (partial order between matrices) 

The partial order < in 9Jl(X) is employed as follows: For Xi,X2 G Tl{X) , > X2 4^ X2 < Xi 4^ 
Xi — X2 > <^ Xi — X2 is a semi-positive Hermitian matrix. 

Definition 3.3 (A : m{D') 971(1))) 

The map A : Tl{D') Tl{D) is defined as follows: First, X E d)l{D'), is regarded as an antisymmetric 
state pi -=^1 j^iyXjj\I){J\ £Ti.-, Then pi is reduced into Ha by the operation p2 :=Trpi £TlA,&nd 
is converted into the matrix representation Y e DJl{D) with the basis D satisfying p2 — J2i jeD ^0 N) I- 
This transformation X t-^ Y is the map A. 



The derivations of this map A arc investigated in section 
Notation 3.4 {€fj) 

For a set X and every i,j € X, (Bfj G 97l(X) is a matrix with entry 1 only at the {i , j) -component and 
elsewhere. For example, for X = {1, 2, 3}, 2 — ^ ^) ■ 



Thus, [efjl/.jGX is 



«£fi 



/lOO 010 OOlx 
/OOO 000 ooo\ 
000 000 000 ' 

000 000 000 
10 10 1 
000 000 000 

000 000 000 

\ooo 000 000/ 

^100 010 001' 



when X is {1, 2, 3} , which wiU be used in this chapter. This example indicates a 3 x 3 block matrix with 
3x3 matrices, thus it is a 9 x 9 matrix. 

Notation 3.5 (A^t)) 

A'-^^-' is defined as a map X ^ K{X^), i.e. a compound transformation for the map A after the matrix 
adjoint operation (Hermitian transpose). Note that A*^^^ operates on Hermitian matrices, as A operates, 
i.e. for a Hermitian matrix X , A^^'>{X) = A{X) because X^ = X . 

In this chapter "map" is a map between matrices. 
Notation 3.6 (identities) 

Let us assume each of 971, 971' is either ofdJlm or m(X). Then id, Id, Id* are denoted as follows: id is the 

mm m,m' m 

identity matrix of 97t, Id is the identity map on 971, Id'^ is the linear map dJl 3 X i-* (Tr X) ■ id G 97t'. 
971,971' will sometimes be dropped, as in id. Id and W^. 



3.3 Propositions and theorems 
Lemma 3.7 

For scalars x and y, eigenvalues of [ Id (^(xA + yA^^^ 

\m{D') \ 

Proof The matrix considered above is equal to S := 



are 



I, JeD' 



(xA + yA(t))(£f;] . For{t,j),{k,l)eD', 

\ J J I, JeD' 



^ («(?,',)(fe,0 ) = 1 ^■)) ((fc' 1 = I ^ ( X ) ^ Kronecker's delta, and A(t) (e^^^-^^,^,) 



1 ( Sji -Su 

2 ' I -S,k S,k 



Observing the whole matrix ^, it is decomposed into the form of the direct sum 
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^ y 



=1 e (2-2 



2-. ^ ^2-2 + |"3 ) wiere 



e 

l<i<j<k<d 



(j,A;)»i 



fo 1 -A 

1 1 
V-i 1 0; 



-2=0 



l<i<d 



(i,i+l)(S)i 



V-1 



1 -1 

1 -1 

• -1 1 

• -1 1 



-1\ 



1 / 





l<i<d 



(i,(i){g>i 



/l 



v 



1/ 



|Si iias eigenvaiues — ?/ and |. (f S2 + f^a) has eigenvalues \y and + \y 
Lemma 3.8 

Let\{x,y)=Tn&yi{\-yl\\yl\^x + \y\} . Then 



arg minA(a;, y) 

x+y=l 



1 rf- 1 



d' d 



and 



min A(.,,) ^ A(i.^-1^ " ^ 



Notation 3.9 (A, A) 

Denote that 



A 

A = 



d-l 



and 



iA + ^A(t). 
d d 



Note that due to the last two lemmas, 



-Aid < I Id (g)A 



< Aid, 



i,.JeD' 



□ 



(3.1) 
(3.2) 



(3.3) 
(3.4) 

(3.5) 



i.e. the absolute values of all eigenvalues of the central side are not larger than ()3.3|) . Note that (|3.4|) 
operates on Hermitian matrices, as A does. 
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Notation 3.10 

Z?'" indicates the bases index ofTiP^ , if D' is the bases index ofH^, as it associates the direct sum of 
the set D'. 



Lemma 3.11 



Id 

OT(D') 



1 f A" Id* -A®") 



> 0. 



i,.JeD' 



Here, Id* is a map from 371(1)')^" to 071(1))®". 
Proof The inequahty H3.b]) is equivalent to 



(id® A" Id*) 



> ( Id«)A«^" 



The following is enough to show the above inequality. 



LHS. = A" 



id# ((sf;-) 



= A"['5/j id],^j-A"id 



RHS. 



( Id^" ^A®") f \<tf'] 
\m(D') JV 



( Id ®A)[ef;] 



< Aid =A"id (< due to ). 



(5n 



(3.6) 



(3.7) 



(3.8) 



□ 



The last lemma successively induces the next two propositions. 

Proposition 3.12 

A" Id* -A®" is a CP map. 

This is due to (ESJ and 1 ^'*^°' ^^ 1. (Any map T : a7t„, OJl„ is CP if and only if 
(Wot™ «)r) [(£ji], ,,^1 ,„ = [r (>£ji)],_j=i...„ is a semi-positive matrix.) 

Proposition 3.13 

A®" (a:) < A" id for X e m{D')'»" if X > , Tr AT = 1. 

Proposition 13 . 1 31 is applied to the proof of Theorem l3.14l 
Theorem 3.14 

E{\'^))> nlog2 ^ for any pure state l^-) G Ti'^". 

This is because Proposition 13 . 1 indicates that all eigenvalues of the reduced matrix from any antisym- 
metric states are less than or equal to (hzt) 

Lemma 3.15 

£'f(c) > 7^-log2 for any density matrix a supported on Ti^" . 
Proof ATote that the entanglement formation is defined as 

Ef{p)= min ^p,E{^i) 
(p.,|<i..>)^eA(p) , 



(3.9) 



where 



Hp) = > 0, ll^.ll = m.Y.P^ = i,^p.|<i>.)(<i>.| = p| 



(3.10) 
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is the collection of all possible decompositions of p. It is known that all of |<I>j) induced from A(p) satisfy 
e Range(p) , where Range(p) is the image space of a matrix p, which is a collection of p\4>) with 
lip) running over the domain of p. Hence 



Ef{p) > min{£;($)|$ G Range(p), ||<I>|| = 1}. (3.11) 

The condition of the lemma above implies Range(p) C iJ^" ^ therefore the last theorem implies Ef{<j) > 
n. □ 



Due to Ecip) = lim — - |Hayden-Horodecki-Terhai oi|^ ^^16 entanglement cost is given as follows: 

n^oo n 

Theorem 3.16 

Ec{o') > ?ilog2 -jTj for any density matrix a supported on 7^*5". 
Corollary 3.17 (The lower bound of the entanglement cost for 7Y ) 

Ecicr) >log^^ (3.12) 
for any density matrix a supported on 7i_ . 

3.4 Conclusion and discussion 

Tliis chapter gave a lower bound of tlic entanglement cost of antisymmetric states for d-dimensional 
antisymmetric states as the inequality H3.12|l . 
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Chapter 4 

Additivity of Entanglement of Formation of Two 
Three-Level Antisymmetric States 



The additivity of Ep is an important problem. The feasible problems to reach a hint of a breakthrough 
were sought. The author tried the problem of 3 (g) 3 plus 3 (8) 3 for antisymmetric states. 

This chapter again focuses on antisymmetric states. The author proves the entanglement of forma- 
tion is additive for any tensor product of two three-dimensional bipartite antisymmetric states. 

This chapter is a derivation from Ishimono 02 |_ 



4.1 Antisymmetric states 

Let us start with an introduction of the notations and concepts used in this chapter. 7i_ denotes 
an antisymmetric Hilbert space, which is a subspace of a bipartite Hilbert space Hab = 'Ha ® "Hb^ 
where both Ha and He are three-dimensional Hilbert spaces, spanned by the vectors TC- is a 

three-dimensional Hilbert space, spanned by states {\i, =23, 31, 12 , where the state \i,j) is defined as 
l^^l^^^-^^l'^ in this chapter. The space 7i_ is called antisymmetric because by swapping the position of 

two particles of states 1-0) in TC-, we get the state —\ip)- 7i?" is the tensor product of n copies of TC-. 

(i) (?) 
These copies will be distinguished by the upper index as H_ , with j = 1, . . . ,n. We assume H_ is an 

(i) (i) 

antisymmetric subspace of '^s ■ 

4.2 The result and proof sketch 

It has been shown in i vidai-Diir-ciri Toil ^Yia,t Ef{p) = 1 for any mixed state p E B{H-). This result 
will play the key role in our proof. Here we will prove the following theorem: 

Theorem 4.1 



Ef{pi ® P2) = Ep[pi) + Ef{p2) 2) (4.1) 



for any pi, p2 G B{'H-). 

Proof To prove this theorem, it is sufRcient to show that 



Ef{pi<^P2)>2 (4.2) 
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since the subadditivity Ep{pi ® P2) < Ep{pi) + Ef{P2) = 2 is trivial. Indeed, it holds 

Ef {pi ® P2) - inf 

< inf ^pfVf^i?(|v^f))(v.f 'I ® 

+ inf^pf)i?(|V;f^)(^f)|) 
= ^f(pi) + Sf(P2) (4.3) 

where subject to the condition pj — X]iPi"'''lV'F^)(V'i"'^l- To prove 1^4. 2}i . we raise the 

following proposition and prove it. 

> 2, for any pure state jV-) e H'^^. (4.4) 
Using the Schmidt decomposition^ , tiie state 1-0) can be decomposed as follows: 

where pi,p2,P3 > 0, = 1, and {|-0|"''')}f^;^ is an orthonormal basis of the Hilbert space h'^^ 

for j — 1,2. JVote tiiat this Schmidt decomposition is with respect to : 'H^^\ or, ^7i^' ^H.^^^ : 

{n^A^ ® H^^) . It is not {n^l^ ® n^l^^ 

■ {^^b' ® ^s'') • Here ":" indicates how to separate the system 
into two subsystems for the decomposition. 



Fact 4.2 

If is an orthonormal basis of 7i„, then there exists a unitary operator U, acting on both Ha 

and Hb, such that U ®U maps the states IV'2), IV's) ^'^^o the states |2, 3), |3, 1), |1, 2), respectively. 

The proof appears in the next section. Because of Fact \4.'A there exist unitary operators J7'^\ [/'■^•' such 
that 

(C/W®C/W®[/(2)®[/(2))|^) 

= E ^^P^\hJ)®\^,J) , (4-6) 

ij=23,3ia2 

where P23 = Pi, P31 = P2, P12 = P3- We denote with the value of 1^4.6]) . 
As is written in the following, we use the following fact: 



Fact 4.3 



i?(|^')(V'|) > 2 if 



P23,P31,Pl2 > 
P23 + P31 + P12 = 1 



(4.7) 
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The proof of this fact also appears in the next section. Local unitary operators do not change the von 
Neumann reduced entropy. Thus E{\ip){ip\) — E{\ip'){tp'\) > 2. Therefore the claim 14. 4|) is proven. 
The entanglement of formation is defined as 

Ep{p)^ inf Vft£;(|V;.)(^.|) (4.8) 
[(p.,iA.)]ieA(rt ^ 

with 

and it is known that all induced from A{p) satisfy l^pi) S Range(p), where Range(p) is is the set of 
plip) with lip) running over the domain of p, called the image space of the matrix p. Hence 

Ef{p) > inf {Eili^M) I IV;) e Range(p), (V^|V) = l} . (4.9) 

Since pi ® p2 G S(H?^) and Range(pi ® pz) C is proven. Therefore H4.1]} has been shown. □ 

4.3 Proofs of lemmas in this chapter 

The proofs of two facts which appeared in the previous section are provided in this section. 
Lemma 4.4 

If C 7i_ is an orthonormal basis, there exists a unitary operator U, acting on both Ha and Ti.B, 

such that U®U maps the states |V'i), IV'2), IV's) into the states |2, 3), |3, 1), |1, 2), respectively. Proof Let 
us start with some notational conventions. In the following, -tl stands for the transpose of a matrix, 
stands for taking the complex conjugate of each element of a matrix, and denotes the transformation 
defined later. 

Let U be represented as ^21 «22 "23 ^ with respect to the basis 1 1) , 1 2) , 1 3) . An operator and its matrix 

( ^11 ^12 ^13 \ 

representation might be different objects, but we identify U with ( «2i "22 "23 j here for convenience. 
Lengthy calculations show that when a 9 x 9-dimensional matrix U ^ U is considered as a map from 
?i_ into 7i_ , it can be represented by the following 3 x 3-dimensional matrix, with respect to the basis 
|2,3),|3,1),|1,2) ; 

^U22U33 - U23W32 U23U3I - U21U33 U21U32 - U22U3l\ 

U32U13 - U33W12 U33UII - U31U13 U31U12 - U32U11 ■ 

\Wi2M23 - U13U22 U13U2I - MIIM23 MIIW22 - U12U2I / 



Then one can show that 



C/«-^[/=(detC/)fo?o 
^ ^ Vo 1 



and multiplying from the right in the equation above, one can obtain = (det U) ■ , since U is 
a unitary matrix, and • W is equal to the identity matrix. 

Since {|V'i)}i=i.2,3 is an orthonormal basis of 7i_, there exists a unitary operator on TL- such that 
|'0i) 1-^ |2, 3), \1lj2) ^ |3, 1), IV'3) |1, 2). Denote 8^ the corresponding matrix with respect to the basis 

{N,i)}ii=23,31,12- 

Let := (det 9^)^ • 9^ It holds that = 9.0.3 Therefore U^(»U^^ U^. The operator 
is the one needed to satisfy the statement of Lemma \4.'A □ 
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Lemma 4.5 



2J 



E{\^P'){^'\)>2 if <^ 



ij=23,31,12 



P23,P31,Pl2 > 
P23 + P31 +P12 = 1 

Proof Let P32 := P23,Pi3 P3i,P2i := Pi2- Tiien it iioids, 



1^') = XI \/p^i«'j>i*>i> 

l<i<j<3 
1 ''^ 

l<«<j<3 
1 ''^ 

l<i#i<3 



wiere |iii2;i3«4) denotes the tensor product \ii) (g) \i2) 1*3) <S> , where G 'H^a ' N2) G ^!4^ 
113) € Ti^'' and |i4) G 7i^' , and the condition 1 < i 7^ j < 3 actuaUy means "1 < i < 3, 1 < j < 3 and 

We are now going to calculate the reduced matrix of \t/y){4>'\, denoted S , and decomposed into the 
direct sum as follows: 



/(2) 



Tr 1^1^') (ip' 

i,j.k.l 



\/PijPki Tr 



\ii;jj){kk;ll\ - j j) {kl; lk\ 



l<i^j<3 
l<k^l<3 



- \ij-3i){kk:ll\ + \ij;ji){kl;lk\^ 
= 7 V VPtjPki Tr (|ii;jj)(fcfc;ZZ| + |ij;ji)(fc;;;fc|) 

l<fc5^/<3 

^ XI VP^fePife I")OjI + 4 X ^'■u I*.?') 

l<fe#'<3 

X /P12+P13 VP13P23 VP12P23 \ X 9 1 9 1 

- VP13P23 P12+P23 Vplipl? ® t(P12) © t(P13) ® t(P23) , 

4 V \/Pl2P23 VP12P13 P13+P23 / 4 4 4 



(4.10) 



wiere ® denotes the direct sum of matrices, and denotes the direct sum of n copies of same matrices. 

We need to get the eigenvalues of S in order to calculate the reduced von Neumann entropy 
= ^Tr(Slog2S) = — J2 ^ log2 A. In this case, fortunately, the whole eigenvalues can 

Aic.v.of E 

be determined explicitly from 



1- cos 61 l-cos(6' + ^) l-cos(6' + ^) P12 P13 P13 P23 P23 



6 



6 



6 



(4.11) 



for a certain < < The eigenvalues above are denoted as {Ai, A2, . . . , Ag}, respectiveij. 

Although A4, . . . , Ag are trivial, Ai, A2, A3 are tiie roots of the cubic polynomial 



16 



16 



(4.12) 
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that is the characteristic polynomial function of the cubic matrix that appeared in H4.1Up . We must 
solve this cubic equation to obtain H4.ll}} . The cubic equation g(A) = is a Cardan's irreducible form,^ 
because S is a density matrix. In such a case, the roots of the cubic equation are 



Itt 47r 

a + f3 cos 6, a + f3cos{9 + —), a + /3 cos{6 + — ) 

o o 

3/32 



(4.13) 



One can easily show that Ai + A2 + A3 = 3a, and A^ + A2 + A3 — 3a + |/?~. If Ai, A2, A3 are equal to the 



roots of the cubic equation + ai A'^ + a2A + 03 0, then Ai + A2 + A3 



-fli and A^ + A2 + A§ = of — 2a2 



hold. Taking ai — — ^, 0.2 — jq from \4.1'J}^ . we get the simultaneous equations {3a — 5, 3a^ + |/3^ = |}, 
with one of the solutions = — Applying this argument to H4.13\} . we complete H4.ll}) . 

The author's idea is now to show that 

9 

E{WM\) = E(-^'' log2 A.) > 2. (4.14) 



i=l 



This will be shown if we prove that the following: 



^{-K log2 K) > 1 and ^(-A, log2 A,) > 1. 



(4.15) 



The second inequality is easy to verify by simple calculations. Therefore, to finish the proof of the lemma 
we need to show that 



^(-Aaog2 AO > 1. 



(4.16) 



i=l 



Without loss of generality, one can assume 9 £ [^t'^]-^ Clearly, Ai G [Oj ti] ^-nd A2,A3 — 
I _ Xi±^Xi 3Xj_ ^ |^_^ I j ^ ^ \^ (.g^jj 5g regarded as the solution of the following simultaneous equations: 
Ai + A2 + A3 = i , Af + A2 + A§ = i . ) One can also show that 



-z\og2 z > 



(log2 12) z 

1 I log 4-1 



4(z- 



1^2 



ifze [0,^] 



(4.17) 



2 ' log^ 2 4^ y 4' 

(see Fig.l ). The first inequality of i)4.J7|) is easily confirmed. One way to prove the second inequality is as 
follows: Let f{z) := (— zlog2 ^) ^ (^5 + '°iog'^2^^ (-^ ~ 3) ~ '^^^ ~ 3)^) ■ Diff('i'(^ntiating this expression with 
respect to z once and twice respectively, we can get the increasing and decreasing table as Table. 14.11 
The table indicates f[z) > for z e [j^, i] . Now we indeed get the lower bounds with polynomial 



z 


1 

12 




1 

8 log, 2 




1 

4 




1 

3 


/(^) 


+ 




+ 


\ 









/'(^) 













+ 




/"(^) 









+ 


+ 


+ 





Table 4.1: The increase/decrease table of (— zlog2 2;) — ^5 + '°iog"^2 ^ ('^ ^ 3) ^ ^("^ ~ 3)^) 



functions. 

Combining all of the above inequalities, we get 1^4.16}) as 

-gA,log2A,>H-(^i^|l±^-2)Ai+4A?>l . 
Therefore i)4.15|) and i)4.14|) are successively shown, and our proof is completed. 



□ 
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0.05 
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0.15 
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0.25 



0.3 



z 



Figure 4.1: — 2log2 z is lower-bounded with the two polynomial functions as 




4.4 Conclusion and discussion 

The additivity of the entanglement of formation for two three-dimensional bipartite antisymmetric 
states has been proven in this paper. 



1 Any pure state |^) of a bipartite space 'Ha ® "Ws can be represented as = ^iVM't) ® 
where lipf) are orthogonal bases of Ha, I'fpf} are orthogonal bases of Bb-, and {pi}i is a probability 
distribution. 

2 In the definition above, it does not matter which of the two roots of det 0^ are taken. 

3 Indeed, = (detZ7^)C/^ = (det 6^)1 det 6^ • ((det ev,)^)tev, = 6^, • Note that detU^ = 
det(det(0^)3 6^) = (det 6^) 5 det 6^ because is a 3 x 3 matrix. 

4 The exact value of 9 has no importance for us in the proof. 

5 A cubic equation is said to be in Cardan's irreducible form if its three roots are real numbers. 

6 The set of {Ai}?^^ does not change if 6 is replaced by — ^. Thus we can change the assumption 



Endnotes 




3 ' 3 



] into e [0, f ] . 
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Chapter 5 

A Gap Between Ep and Er, 



Are there gaps between Ed and i?c? That was an important problem. The answer was yes. We found 
such examples for a 2 (g) 4-level system in 2004, after the bound entanglement for a 3 (8) 3-level system 
was found j Horodecki 98 | Poland. 

This chapter is a derivation from iMatsHmoto-ghimono-Winter q S. 

5.1 Background 

The problem is stated as follows: 

Is there a gap between the two important and not-yet-fully-analyzed quantities of quantum 
entanglement presented below? 

1. Ec{p), the entanglement cost, 

— the asymptotic quantity of Bell states needed to produce p under LOCC operation, 

2. Ed{p)^ the entanglement distillation, 

— the asymptotic quantity of Bell states retrieved from p under LOCC operation. 

If there is no difference, meaning Ec{p) — Ed{p) for any p, the problem of quantification of quantum 
entanglement may end up easy to handle. That the equality holds for pure states is already known 

jPopescu-Rohrlich 97|IVidal OOllNielsen OOI|Donald-Horodecki-Rudolph 02| |^}^^g problem resideS for p being nOn- 

pure state. If there is a difference, which means Ec{p) > Ed{pY , one can say the irreversibility occurred 
under LOCC operation and that the degradation of quantum entanglement is inevitable. 

Analyzing this problem is theoretically interesting when investigating the property of quantum 
entanglement. 

Here is a historical note. There are already answers in that there exists such p that satisfy En (p) < 
Ec{p)- The first examples are such states p that satisfy Ec{p) > and Ed{p) — 0, which are named 
bound entanglement j Horodecki ^ _ r^j^^ bound entanglement appears in a 3 (g) 3-level system. Other 
attempts are performed at jV'dai-Diir-girac qi| antisymmetric states. This chapter gives the examples 
of a 2 ® 4-level system, with mathematical software Maple, through calculating the Holevo capacity of 
unital channels 1^'"*^ "^^=1. 
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5.2 Introduction 



We need to estimate both the values of Ec and Ejj . Ec can be calculated by the Holevo capacity if 
the full capacity C is determined through Stinespring's relation because of the relation between Ep and 
the Holevo capacity (see Sec. I1.5|l . To evaluate the value of Ejj there is a useful estimation theorem. 

Notation 5.1 (the logarithmic negativity) 

The logarithmic negativity E^ of p where p is shared by two sites is defined as follows: 

En{p) = \ogY,\K\ (5.1) 

i 

where {Xi} are the eigenvalues of where is the partial transpose of p. The partial transpose of a 
state for a bipartite system is dehned as {ai(g,j^k(g,i){i(»j),{k(»i) (afc<8.j,i«.i)(i<»i),(fc(»0 • 

Theorem 5.2 (tSE^WEEEEM) 

Ed{p) < En{p) (5.2) 

Note that Em is simple to calculate even though the method to calculate Ed for general cases is 
not known. 

Therefore, the aforementioned gap Ed{p) < Ec{p) is concluded, if 

- Ec{p) is calculated and 

- Efq{p) < Ec{p) is determined, 
by some method. 

5.3 Unital channel and its associated state 

We considered a generalized depolarizing channeP on a qubit: 

s—0,x,y,z 




Distillation 



Figure 5.1: Does the irreversibility occurs during the production and the distillation? 
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with the Pauh matrices 




and a probabihty distribution {ps}s=Q,x.y.z- 

There is an interesting theorem about the additivity of the Holevo capacity for depolarizing channeh 



Theorem 5.3 (jKing 02a|) 

For generalized depolarizing channels, additivity of the Holevo capacity in terms of tensor product with 
an arbitrary channel holds. Namely, 



C(A(g) A') = C(A) + C(A') 
for any generalized depolarizing channel A and any channel A' . 



(5.4) 



Note that up to unitary transformations on the input and the output systems, each unital qubit 

channel has this form iKing-Ruskai Ol||Fujiwara-Algoe^^ 



We assume the condition 

PO + Pz - Px - Py > \P0 + Py - Px ~ Pz\,\PO + Px - Py - Pz\, 

which does not lose generality. Then for such a channel A the capacity is given by 

C(A) = 1 - Snuni^) 



(5.5) 



(5.6) 



where S'min(A) is the minimum output entropy of A, which is achieved at either of the inputs |0), 
thus S'„iin(A) — S'(A(|0)(0|)) — S'(A(|1)(1|)) . The optimal ensemble of input signals of the channel in 
order to achieve the Holevo capacity is the uniform distribution (1/2, 1/2) on the two states |0), |1). 

A Stinespring dilation for this map A : B{TL±) — > B(TLo) is given by an isometric embedding U : 
Tii — > Tlo (8) Ti-a where 7ii,7Yo = 'C'^ and Ti^ — C^, in the following 8x4 block form: 



U 



/p^CTx 

fP^Oy 

\ VP^CTz ) 

and the corresponding subspace K C Tio ® '^a is spanned by 



(5.7) 



IV-a) = ^/P^|0) ® |0) + VP^|1) ® \x) 



rpM®\z), 

fp'z\^)®\z). 



(5.8) 
(5.9) 



(You can easily verify that A is equal to the operation p ^ Tr-^^ UpW .) The average of the optimal 
input state"^ is transformed to the equal mixture of the two pure states 



PA^^IV'AXV'Al + ^IV^fX^fl 



(5.10) 



by U. 
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With Theorem [TTTl (|^ leads to 

E^--^{pa) ^ S^in{A) ^ H{po + p„p^ + Py) (5.11) 

where H is the entropy function and Ep^{p\ (g) cr) = Ep^{p/C) + Ep^{(j) for any a G Oa', with an 
arbitrary channel A'. Theorem OS leads to Ef-^'ipf") = nEf-^ip), thus 

ES--^{pA)^Hipo+p,,p.+Py). (5.12) 

This leads to that the decomposition of p'^" into the 2" equally- weighted tensor products of |i/'a)(?/'a|, 
IV'a)(V'aI being is formation-optimal^ . By the convex roof property of Ep this implies that any con- 
vex combination of these states is a formation-optimal decomposition. (This argument was also used 
in [Voiibrecht-'Werner qJ extcud the domain of states with known entanglement of formation.) In partic- 
ular, we can conclude that 

ES'-^ip) = EPip) = H{po+p,,p,+Py) (5.13) 

for any mixture p of |?Aa)('0a| and 

— We may omit the superscription ( ) from here. 



5.4 A gap between Ec and Ej^ 

One can verify that the partial transpose p\ of the optimal state pA is decomposed into a form of 
the direct sum of two 4 x 4-matrices from (|5.10|) . which have the same characteristic equation 



f(2t) = 0, with 

f{t) -t^ + AipaPxPy +PoPxPz +P0PyPz +PxPyPz)t ~ IGpoPxPyPz 



(5.14) 
(5.15) 



This equation /(2i) = has only one negative root denoted Iq, and / is decreasing in a neighborhood 
of io- See Fig.lO 



1 0.04 - 

toi 




V 0.02 - 

, , 1 . „ 




1 -0.8 -0.6 -0.4 -0.2 


, 0.2 0.4 0.6 0i8 1 

/ V 


\ -0.02 -i 





Figure 5.2: The graph of f{t) for {po,px,Py,Pz) = (l^i ll) as an example 



Assume p^ has eigenvalues {to, to, ti,ti, i2, ^2, ^3, is}- Since to -(- + + + ^2 + ^2 + ^3 + ^3 = 1 , 

IIPaIIi = \to\ + \tQ\ + \tl\ + \ti\ + \t2\ + \t2\ + \t3\ + \t3\ = \2to\ + \2ti + 2t2 + 2t3\ ^ (-2to) + (1 - 2to) - l-4to, 

where ||-||i is one-norm, or the sum of the abstract value of the eigenvalues. Hence i?Ar (pa) =log(l — 4to). 

Thus Ejm{p\) < Ec{pa) is successively equivalent to log(l — 4to) < Et{po + Pz,Px + Py), 1 — 4io < 
2mPo+P.,P.+Py)^ i-2«(^o+P.,P.+P.) ^ 2i„, and 



OH{po+Pz,P:,+Py) _ 1 ~ 

f{-- . ^1 >o. 



(5.16) 
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V \^ ( E d(p) < E c(p) 
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Figure 5.3: Pursuit oiEu{p) < Ec{p)'- The issue resides whether there exist points satisfying H5.16|l with 
(|5.5|) . The resuh is, at least aU the points between the front of the curved face inside the pentahedron 
are such points. Here, the curved face corresponds to H5.16|l and the pentahedron corresponds to H5.5|l . 



This means if Px, Py, Pz satisfy this incquahty H5.16|l . there is a gap between the entanglement 
cost of PA, and its entanglement distillation, namely E^Ipa) < Ec{ph)\ Fig. 15.31 shows the plot of the 
region of these {p.j;,Py,Pz)- Due to this continuity, also for a mixture of \iPa){'4'a\ and |^a)(V''a I which is 
sufficiently close to p\, we observe a similar gap. 

Example 5.4 

Assume {po,Px,Py,Pz) = (5, ^, ^, Then Ed{pa) < Ec{pa). 

Proof A short calculation reveals that ||pa||-^ = 5/3, so E]j{p\) < log(5/3) w 0.737, which is smaller 
than the entanglement cost Ec{pi) = i?(l/3, 2/3) w 0.918. 

You can also verify that LHS of 1^5.16}} is 0.00784 (which is larger than 0), thus E^Ipa) < Ec{pa)- 

□ 

In this case, 
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(5.17) 
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Example 5.5 

Assume {po,p^,py,p^) = (|, |, ^) with u,v e [0, |)U (i, 1]. Then Ed{pa,s) < Ec{pA,s) holds for 
PA,s = s|V'a)(V'a| + (1 ~ •s)I^a)(V'a I with < s < 1. In this case, pA,s is an 8 x 8 matrix with variables 
s,u,v, represented as 



( su 



sV uu' 





s'v 









s'v' 







su 



s'^/uv s'y/~uv' 

sy/—uv' s^J—u'v' 

Sy/uv s\fu/v 

\ -s'V^ -sV-w'"' 







\ 



—s'\/—uv' 




-s\/—uv' s^/uv 



sV-w'«' 

-s\/—u'v' sVu'v 



S It 




-s'Vwj/ 









S\/—VV' 



■'S\/—vv' sv 









(5.18) 



where s' — 1 — s,u' — 1 — u,v' — 1 — v. 
Proof 

By eq. 1^5.13^) . Ec{pa,s) = 1- The key observation is whether logUpJ^lli < EdpA.s)- This is true 

because the condition ^5.161) is always satisfied, as LHS of H5.1b]) is + i)^ — (u— — i)^ 



which is zero at {u, v) — it is also an increasing function w.r.t. both u = (m — and v = (v 

foru,ij < i (SeeFig.\^. 



1\2 



□ 
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0.12^ 
0.08 
OM 

04 




Figure 5.4: The graph of i(u - i)^ + l{v - i)^ - (u - i)2(„ _ i)2 



5.5 Conclusion and discussion 

We have found the states which show the gap as Ed {p) < Ec (p) for a 2 4-level system. They 
were found by analyzing general depolarizing channels with Stinespring's dilation. This discovery also 
shows a usefulness of the relationship between Ep and the Holevo capacity (see Section [T3|) found by 
jMatsumoto-Shimono- Winter 04 | ^ rpj^^ examples showu in this chapter are the lowest dimensional cases among 
those bipartite states p being proved Ei:,{p) < Ec{p), in that they are 8 dimensional (2 (g) 4), though 
the next lowest cases are the bound entanglement which are 9 dimensional (3 ® 3) found in earlier in 
jHorodecki ^ ^ Q^j, research may exemplify that analyzing the entanglement even in the easiest cases 
would involve variant ideas such as quantum channels that are tied to entanglement by Stinespring's 
relation. 

Endnotes 

1 Ec{p) < E£i{p) docs not occur for any p in principle. If it did, quantum entanglement would 
increase infinitely by LOCC operations. 
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2 A depolarizing channel and a generalized depolarizing channel: 

- A depolarizing channel for a qubit is a channel with its output ellipsoid being a sphere and its 
center located at the center of the Bloch sphere. ( The output ellipsoid is explained in Subsection 



- A generalized depolarizing channel here is defined as follows: a qubit channel is a generalized 
depolarizing channel as long as the center of its output ellipsoid is located at the center of the 
Bloch sphere. 

The output ellipsoid of the generalized depolarizing channel of H5.3|l has three axes in x,y,z directions 
with the radius of Po + Px - Py - Pz,Po + Py - Px - Pz and po+ Pz - Px - Py 




Figure 5.5: An illustration of the Bloch sphere and the output ellipsoid of a general depo- 
larizing channel. With the condition H5.5I) . the longest axis is z-directional. 

3 Consider the definition of the Holevo capacity: 



We call the optimized {pi\i the optimal input. The average input is '^^PiPi- 

4 The formation-optimal here means that the collection of states gives the optimal value of 11.24|l , 
which is the formula of the entanglement of formation. 



m ) 



n 



n 



C{A) 



max 

n-pi,...,p„;pi,...,p., 
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Chapter 6 



Numerical Verification of Superadditivity 



The strong superadditivity of Ep is rather simple to check among the equivalent additivity problems. 
The author numerically checked the lowest-dimensional cases, which are described as "2 (g) 2 plus 2 2". 

In this chapter, the strong superadditivity of the entanglement of formation is numerically verified. 
This is a part of attempts to verify the expects that the whole additivity holds for the entanglement 
of formation (see Sec ll.6|) . It is interesting because if there were any single counterexample against the 
strong superadditivity, the expected relation Ep = Ec and the additivity of the Holevo capacities would 
collapse. The additivity of the minimum entropy of channels would also collapse ISES^Hl, 

As a result, any single counterexample has been found from approximately a million randomly- 
chosen cases and various directive searches. Even though this is a negative result, it gave the knowledge 
that finding counterexamples is a challenging problem even if they exist. 

6.1 Introduction 

The strong superadditivity of the entanglement of formation for four-partite qubits is numerically 
checked in this chapter. The strong superadditivity is represented by the inequality defined in Section 
1121 as 

Ef{p) > Ef {Til p) + SF(Tr2 p) 

for p on Hai^'Ha2®'Hbi®'Hb2- (6.1) 

The additivity of entanglement of formation can be deduced directly from this superadditivity if 
this inequality JHIIl holds for any pure case iMatsumoto-Shimono-v^inter ^. 




Figure 6.1: The left hand side of 16.11 corresponds to left and the right hand side corresponds to right. 
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In this chapter 7iAij^yi2i^si and 7is2 a-re assumed to be a two-level systems (four qubits). There 
are several reasons. First of all, this is the lowest-dimensional case. Second, there are feasible method 
to numerically calculate the values appeared in (|6.1() : 

• The left hand side of (|6.1I) is just a von Neumann entropy when p is a pure state. 

• Each term of the right hand side of (|6.1|l is feasible to calculate by the method of the concurrence 
jmii-wootters 9 7|| Wootters 98 ) ^j^g^^ argument is as low as two-dimensional. 

Note that direct calculation of Ep is formidable as it is defined as an optimization form on high- 
dimensional space as it is defined as p.24|l . 

Furthermore, there is another reason to check the inequality rather than to check the additivity. 
Technically, inequality check is easier than equality check. You cannot stop concerning about the cal- 
culation accuracy when you check the equality numerically, whereas you can when you are checking the 
inequality numerically unless the equality exactly holds. 

6.2 Calculation procedure for the inequality 

In this section, we show the procedure how to numerically verify H6.1|l for a chosen case. This is 
the base of the experiment executed millions times, which works for Ti^ij , Tis^ , Tiy^^ ; ^-82 being 2-level 
systems. 



(1) Start from a four-partite qubit pure state that is represented as sixteen-dimensional vector 



= (o^oooo ooooi • ■ • aiiio aim 



(6.2) 



for a randomly chosen 16-tupIe of complex numbers {aijki} with the normalized restriction 



i 



(2) Calculate three reduced density operators, PA1A2, PAiBi and PA2B2 from 
I^AiSiAsSaX^AiSiAaSa I , as 



PA1A2 = TrSjS2 |*AiSiA2S2X*AiSiA2-B2i, 
PAiBi = TrA2S2 j*AiSiA2S2X*AiSiA2-B2j, 
PA2S2 = TrAiSj |*AiBiA2-B2X*AiBiA2S2l- 



(6.3) 
(6.4) 
(6.5) 



12 3 
4 5 6 7 
8 9 A B 
_C D E F _ 

is the summation of the elements of pa^b^ = Tiai^s |^AiBiA2-B2)(^Ai_BiA2-B2 I positioned by the 



They are 4x4 matrices. Assume they are represented as 



Then each symbol 0, 1, . . . , F 



same four letters as follows. H6.3^ . H6.4fl . H6.5fi for 



2 . 3 

.0.1 2.3.... 

4.5 6.7 

.4.5 6.7.... 

....0.1 2.3. 

0.1 2.3 

....4.5 6.7. 

4.5 6.7 

8.9 A . B 

.8.9 A . B . . . . 

C . D E . F 

. C . D E . F . . . . 

...-8.9 A - B . 

8.9 A . B 

. . . . C . D E - F . 

C . D E - F 



12 3 

4 5 6 7 

& 9 A B 

C D E F 



12 3 

4 5 6 7 

8 9 A B 

C D E F 



12 3 

4 5 6 7 

8 9 A B 

C D E F 



respectively. 



12 3 

4 5 6 7 

8 9 A B 

C D E F 
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(3) Calculate the Ef for a pure state A1B1A2B2) , as 

£F(l*AiSiA2-B2)(*^li3lA2S2|) =S{pAiA2), (6.6) 

where S{paiA2 ) is the von Neumann entropy that is calculated by the Shannon entropy — Ai log Xi 
with {Xi}i the eigenvalues of the reduced density operator PA1A2 that is a 4 x 4 matrix. Thus one 
need to solve a quartic equation (a polynomial of order four) or utilize some mathematical package 
to calculate the eigenvalues of matrices. 



i 



(4) Calculate the 'spin-flip' transformations | Wootters _98j ^j- ^j^g reduced density operators, PAiBi 
and PA2B2 by calculating 

PA^B^ = (o-y®o-i/)PAiSi('7i/®cry), (6.7) 

PA2B2 = {(Jy(^(jy)p^A2B2{'^y'^<^y)- (6.8) 

/ -1 \ 

^ ) is one of the Pauli matrices. Thus ay ® Oy — [ n ? n n I • This spin-flip 



/ illUB w - 10 

' V -1 

F -E -D C 



transformation p 1-^ p is the operation like 



12 3 

4 5 6 7 

8 9 ,4 B 

C D E F 



-B A 9-8 
-7 6 5-4 
3-2-1 
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(5) Calculate the concurrence of pa^Bi and PA2B2 by the following calculations: first calculate the 
square roots of the eigenvalues of pAiBipAiBi and arrange them in decreasing order Ai, A2, A3, A4; 
next, calculate the concurrence 

C(pAiSi) = max{0, Ai - A2 - A3 - A4}. (6.9) 

Do the same for PA2B2, as 

C{pA2B2) = max{0, pi — P2 — P3 — Pi} (6.10) 
where pi,p2,P3,P4 are the eigenvalues of PA2B2PA2B2 arranged in decreasing order. 



i 



(6) Calculate Ef for pAiSi and PA2B2 by 



EEiPA.B.) = I ^~^^~':^p-^-^y- 1 , (6.11) 



Ef{pa2B2) = I ^-^^-CiPMB2y- I ^ 



where H2 is the binary entropy function. 



^Ef{p) = ^ ^-^V^ iog, _ i±%AElw!iog, i+,A7^ ) 



(7) Compare the values 

£F(|*AiSiA2B2X*-4iSiA2-B2i) (pAi Si ) + (PA2S2 ) , (6.13) 

to check whether the superadditivity of Ef holds or not. 
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6.3 Procedure to pick up random points and their results 

Three different metlrods are performed to pick up points, wliicli are tlien substituted into the pro- 
cedure described in the section above. 

6.3.1 Random search 

Sixteen parameters Ziji^i (i,j,k,l = 0,1) are chosen from uniformly from the square whose four 
vertices are ±1 ± \/— T on the Gaussian plane. Then they arc normalized by Uijki = 



These values are substituted into A1B1A2B2) = X); j fc ; Q^ij'-fe-d*) li) 1^) 10- 
Fig. 16.21 shows a typical case. Each of the dots is located at 

{x,y) = {S{paiA2),Ef{pAiBi) + EpipA^B^)), 

The result is that there has not been any counterexample for approximately million points. 
2 I r 



(6.14) 



1.5 - 




0.5 - 



Figure 6.2: Randomly chosen 10,000 states from the whole four-qubits space. There is not any exception against 
the inequality, as all the points in the figure satisfy the condition x < y. 



6.3.2 Pseudo zero-neighborhood search 

If \'^AiBiA2B2) is a separable state with respect to one side AiBi and the other side A2B2, equality 
holds for the inequality. Thus to find the counterexamples, it would be worth to check the neighborhood 
of these separable states. We employed the following method to pick points around the separable states. 
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Set 

i*> = E"yi*)ij')®E"'='i^)io+e J2 a..,,kAm\m) (6.15) 




Part of the results are shown in Fig l6.3l for e — 0.2 and 0.05 with 1,000 points for each. The author has 
checked for many parameters e, but there have not been any of the aforementioned counterexamples. 



Figure 6.3: 1,000 points are randomly chosen from each of certain neighborhoods of the separable states 
for each figure, left: e — 0.2, right: e — 0.05. There is also no exception against the inequality, as all 
the points in the figures satisfy the condition x < y. 

6.3.3 Minimum search 

Seeking the minimum value of 



is also performed. The issue is whether the value of ()6.17|l can reach below zero. The employed procedure 
is as follows: 

• First, pick a random point. 

• Next, repeat the following step: 

— pick hundreds of points in a certain neighborhood of the previous points and 

— the point giving the minimum value of (|6.17|) is chosen to be the next point, 

until the movement of the points is negligible. The neighborhood shrinks as the steps progress. 

The result is that there has not been any counterexample for the inequality for tens of trials (see 
Fig. 16.41) . We observed that all of the minimum values of these experiments are zero with the convergence 
in the separable state represented by \'i'AiBi){'^AiBi \ ^ I^A2B2)(^A2B2 I- 




EpipAiBi) + Ef[pa^B-2) - S{pAiA2 



(6.17) 
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o' ^ ' ^ ' ^ ' ' 

0.5 1 1.5 2 

Figure 6.4: The trails of seeking the minimum. Seeking the minimum oi x ~ y was tried tens of times. 
Each trial is distinguished by color. Every local minimum was zero. 

6.4 Results and Discussion 

6.4.1 The experiments and the result 

More than one million cases for the random search and tens of the directive search described above 
were performed. As a result, no counterexample was been found. 

6.4.2 Review of the search ways 

The experiments described in this chapter were a pilot test and there are many points to be improved. 

For the random search, sixteen parameters {zijki} would prefer to be chosen from the Gaussian 
distribution rather than squares of the Gaussian plain, to remove the bias on the pure states. 

At first glance, one million points seems almost myriad. The space those states reside on is, however, 
not two- or three- dimensional space. The space \'^AiBiA2B2) resides on is essentially fifteen dimensional 
in real^ . Since ylOOOOOO w 2.5, a figure with the diameter of a third (^) of the whole space would not 
be detected by a million random detections, in a rough estimation. Thus a numerical check performed 
here must be treated with care; if the region of the counterexamples has a diameter less than a third of 
the collection of the whole pure states, such a region is not detected. 

For the (pseudo) zero-neighborhood search, we observed that on the plots (see Fig l6.3(l . the points' 
distribution is strongly determined by the parameter e of H6.15|l , as the possible locations of the points 
on appearance do not overlap. This fact may bring us some analytical knowledge. 

It would be better to employ a more refined way to seek the minimum values, as there are many 
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methods for seeking minimum values, though we saved labor to implement the program. Then we can 
run a lot of trials to seek the minimum value. Since the result to seek the minimum value of (|6.17|) shows 
that every local minimum is zero for tens of experiments, even if the minimum were less than zero, the 
minimum is hard to reach. 

The experiments described above might some hints for the future analytical proof about whether 
the inequality holds or not. 

6.4.3 For higher dimension 

The experiment performed above is for the lowest dimension that is two for each of the four systems. 
This results in turns to be the substituted pure state being sixteen dimension (2"*). What would happen 
in higher dimension? Even for the next lowest dimensional case, we do not know a feasible way to 
calculate the entanglement of the formation, contrary to the case of two dimension. That makes the 
problem difficult. 

Endnotes 

1 The pure states form a thirty-one dimensional manifold. Unitary operations, which is four 
dimensional, on each four of qubits does not effect any difference on the inequality. Thus 31 — 4x4= 15 
is the virtual dimension. 
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Appendix: Maple 6 program 



The following is one of the programs of Maple 6 used for the random search (Subsection. I6.8.1|l . 

vn\h{LwearAlgebra ) 
gauss := proc() 
localv7, v2,fac, r; 
r:= 3.0; 

while 1 < r or r = do 

vl :=rand( )x.2*l(r(-ll) - 1.0; 
v2:=mnd( )x.2*10^(-ll) - 1.0; 
r := vlxvl + v2xv2 
end do; 

fac := sqrt(-2.0xlog(r)//'); 
return vlxfac + Ixv2xfac 
end proc 

generator := proc() 

Vector(16, eval(/^ gauss( ))); %'sqrt((HermitianTranspose(%)) . %) 
end proc 

/:=Vector([0, 1,4,5]); 

h:= 'if (0 < A', -X log2(A'), 0); 

h2:= hix) + h( 1 - x); 

a := <<0, 0, 0, -1 ) I <0, 0, 1, 0), (0, 1, 0, 0), (-1, 0, 0, 0)); 



judge := proc(a) 

local p, pO, pi, p2, SO, SI, S2; 

p:= a. (Hem]itianTmnspose(a)); 
pO := MatnK(shape = hermitian, 4, 

add(p[/[/]x2 +/[*]+ l,/[/-]x2+/[/t]+ l],yt= 1 .. 4)); 
pi := Matrix (5/; ape = hermitian, 4, 

{i,j)^ add(p[/x4-A:,7x4-A:],A:=0.. 3)); 
p2 := Matrix (5/; ope = hermitian, 4, 

ii,j)^ add(p[i + fcx4,7 + /:x4],A:=0.. 3)); 
Egenvalues(pO); 
50:= add(h(°/3ik]),k= 1 .. 4); 

sort(map(sqrt@5fl, evalf(Eigenvalues(pl . a . (Transpose(pl)) . a)))); 

51 := c2e(imx(0, %[4] - %[3] - ?^2] - ?^ 1 ])); 
sort(map(sqrt@9J, evalf(Egenvalues(p2 . a . (Transpose(p2)) . a)))); 

52 := c2e(rmx(0, %[4] - %[3 ] - %[2] - %[ 1 ])); 

if .65 < (SI + S2y SO thenpikit(SO, SI, S2);prkit(seq(a[i],i= I.. 16)) 
end if; 

printf("(%f,%f)," , SO, SI + S2) 
end proc 

for/ to 250 do judge{generator( )) end do 



The output of this program may be as follows: 

(1.423552, .077738) , (1 . 241545, . 057057) , (1 . 227051 ,. 029807) , (1 . 301131 ,. 075570) , (1 . 366521 ,. 004438) , 

(1.235635, .175308) , (1 . 446036, . 026785) , (1 . 299543 ,. 082059) , (1 .480392 ,. 178842) , (1 .010854, .224507) , 

(1.515526, .209590) , (1 .424525, . 195779) , (1 .478113 ,. 011087) , (1 . 161710 ,. 040721) , (1 .416082, .311309) , 

(1.120448, .072590) , (1 . 528262, . 000081) , ... 
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Part III 
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Chapter 7 



Calculating the Holevo Capacity 



The Holevo capacity, the classical information communication efficiency of a quantum channel, was first 
formalized as the Holevo quantity in f 973 in Russia. The calculation method of the Holevo capacity of 
an arbitrary qubit channel with the guarantee of computation convergence was performed for the first 
time in 2004, which is presented in this chapter. 

This chapter is a derivation of i^^y^^^^~^™^^~^^^^^^°^Q~^^^^^^~^^^^°^° Q^l a,nd 

7.1 Significance in calculating the Holevo capacity 

7.1.1 The Definition of Holevo capacity 

Here we review concepts of the Holevo capacity of quantum channels. 

A quantum channel is a CPTP mapping from its input space 6(7^ in) to its output space S(7iout)- 
(One may supplement the explanation, as this quantum channel is memoryless.) This chapter mainly 
deals with qubit channels, thus both their input and output spaces are 2-dimensional and each channel 
is characterized by its output ellipsoid presented as Fig. 17.21 Upper Left. 

The Holevo capacity for a quantum channel A : B{Ti.in) B{Hout) can be represented in two ways: 

1. One way is jHolevo 73IIHolevo 771 ISchumacher- Westmoreland 971 

n n 

C(A) = max S'( Vp; Ap,) - Vp, S'(Ap,), 

ri; pi,...,p„;pi,...,p„ — ' — ' 

.= 1 ^ .=1 ^^^^^ 

subject to n e N,pi > 0,^pi = 1, pi € B{Hin) (i = 1, . . . 

i=l 

where S{-) is the von Neumann entropy. We call each element pi of the optimized (pi) an 
engaging input of A, and the number n the engaging number. The engaging inputs are 
the input states which are used to send messages through the quantum channel when the 
communicating efficiency achieves the Holevo capacity. 

Note that H7.1|l is the maximization form over the average entropy of the outputs minus the 
entropy of the average of the outputs. Therefore, it is interpreted as, the maximum vertical 
distance between the following two surfaces: i) the epigraph of the von Neumann entropy 
function of which the domain is constrained to the output image of the given quantum channel, 
and 
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Given a function, its epigraph is the set of points on and above the 
surface of its graph. The convex roof is the function determined by 
the convex huU of the epigraph. The vertical distances are indicated 
by the arrows in the left figure. 



Figure 7.1: The vertical distances between the epigraph and the convex roof 





Upper Left : 

The output ellipsoid inside the Bloch sphere 
representing a qubit channel of 

A4 : p{x, y, z) t-* p{0.6x + 0.021, 0.601y, 0.5z + 0.495). 

Lower Left : xz-section of the Lower Left figure 
Green circle: the Bloch sphere 

the minimum enclosing sphere 
w.r.t. the quantum divergence 
the output ellipsoid 
the center of the red "circle" 
w.r.t. the quantum divergence 
In the Lower Left figure, the Bloch sphere (green) and 
the output ellipsoid (yellow) are a sphere and an ellipsoid, 
respectively, w.r.t. the Euclidean metric. The enclosing 
sphere is, however, like an ellipsoid and is subtly distorted, 
even though it is a "sphere" w.r.t. the quantum divergence. 

The color of the surface of the output ellipsoid in 
the Upper Left figure corresponds to the quantum 
divergence H {p' \ \a-') with the following color: 



Red oval: 

Yellow ellipse: 
Red cross: 



0,3D4 0.306 



0,316 0,318 



Here tr' is the center of the minimum enclosing sphere 
(black cross of the Lower Left figure) and p' is substi- 
tuted with points on the output ellipsoid. Therefore, the 
most reddish points are the engaging outputs. 



Figure 7.2: The Holevo capacity as the quantum divergence radius of the minimum enclosing sphere 
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ii) the graph of the convex roof of this function. See Fig. 17. II for a conceptual illustration. 



2. The other way is 



johya-Petz-Watanabe 97| 



C(A) 



min max H(Ap 1 1 Ac 

creB(Ki„) pG8(Ki„) 



(7.2) 



where is the quantum divergence. It is the radius of the minimum enclosing sphere 

containing the output ellipsoid when the radius is measured by the quantum divergence. (See 
Prop. lLBl to refer to the pseudo distance properties of the quantum divergence.) 



The physical meaning of the Holevo capacity is as follows: it is a classical information capacity^ of 
a given quantum channel at which input particles arc not allowed entangled with each other, and the 
output particles are measured collectively. 

To fully utilize quantum aspects of the quantum channel, one may consider the capacity at which 
the input particles are freely entangled with each other to send messages in order to take advantage of 
the communication efficiency. This capacity is called the full capacity, and it is represented as 



7.1.2 Why do we calculate the Holevo capacity? 

For historical background, see Subsection 12. 1.31 

What if the Holevo capacity C(A) can be calculated? The benefits are as follows: 

• The actual benefit would occur when the full capacity C{A) = lim^v^oo C{A'^^)/N is calculated 
for arbitrary channels. (The full capacity provides the theoretical upper-bound of the channel 
capacity in the age of quantum information commencing tens of years from now, when each light 
particle, i.e. photon, carries a significant amount of information.) Therefore, it is very desirable 
that C(A«'^) for = 1, 2, . . . , oo can be calculated. 

• As a matter of theoretical contribution, theorists like to verify whether the additivity holds for the 
Holevo capacity, that is, whether C(Ai O A2) = C'(Ai) + C(A2) holds for arbitrary channels Ai 
and A2. (The author is interested in finding circumstantial evidence as to whether this additivity 
holds or not. A piece of evidence which numerically suggests the additivity appears in Subsection 




C(A»^) 



(7.3) 



N 



ES) 



7.2 Algorithm to calculate the Holevo capacity 



7.2.1 Difficulty in optimization 



The objective function 



n 



n 




(7.4) 



to be maximized has the following two properties: 



1. Convex w.r.t. {pi)i, the set of input states. 
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a concave function with a convex search space 



a convex function with a convex search space 



Figure 7.3: Right: Calculating the maximum of a concave function with a convex domain is feasible because 
the solution is obtained by simple gradual ascending or hill climbing. Left: Calculating the maximum of a 
convex function is, however, rather difficult. The maximum cannot be generally reached by gradual ascending 
because there could be multiple local maximums, which appear at the boundary of the search space. 

Due to the joint convexity johya-Petz 93|^ ^.j^j^^ jg^ 

XH{p\\a) + {l-X)Hip'\\a')>H{\p+{l-\)p' \\ \a + {1 - X)ct' ) (7.5) 

for < A < 1, the quantum divergence H{p\\a) is convex w.r.t. p. Note that the 
objective function 17.4|l is rewritten as 

J2p^H{Ap.\\J2PJ^PJ)■ (7.6) 

Since H{p\\a) (= Ti-{p\ogp — ploga) ) = —S{p) — Trplogcr, as this first term 
is convex and the second term is linear w.r.t. p, the objective function is convex 
w.r.t. ipi)i. 

- 2. Concave w.r.t. {pi)i, the probability distribution. 

The first term >5'(^"^-^ pi Api) is concave w.r.t. (j>i)i due to the convexity of the von 
Neumann function S{-) itself. The second term Y^^i=iPi ^i^Pi) is linear w.r.t. {pi)i. 

In such a function H7.3|l . it is hard to identify the function's global maximum. When one employs a 
method of gradual ascending or hill climbing, there is no way to guarantee whether the obtained solution 
is the global maximum (see Fig. 17.3(1 . 

Because of the convexity of the objective function ((7.3|l w.r.t. {pi)i, each pi of the solution of {pi)i 
appears at the boundary of B{Hin)- Hereafter, we denote the bound of B{Hin) as dB{Hin}- 

7.2.2 Solution — Computation w^ith lattices 

The issue is ascribable to the difficulty in guaranteeing the calculation error. In order to guarantee 
the convergence in calculation, the following method is effective, as it simplifies the problem of maximizing 
a complicated function into one maximizing a concave function on a convex search space. 
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Solve the following approximated problem: 

n 71 

= maji S(S^ Pi Api) -S^ Pi S{Api), 

subject to being a probability distribution, (7-7) 

where £ — {pi}i^i is a set of fixed points covering dB{Hout) 
with a specified (smaU) coarseness S. 



Here the covering dB{TLout) is a set of points, for example, the lattice shown in Fig. 17.41 Left. The 
coarseness 5 is defined as the maximum distance from arbitrary points on the sphere to its nearest 
vertices among {pi}^^i- 

We employ the following lattice, and we denote it ik for a positive integer k: 



The lattice is defined as, 

{pi]^^i being the set of vertices of the lattice which divides both the equator and 
the meridian of the Bloch sphere, as it is a set of 

UTT 2V7T UTT 2vTr UTT \^ 

sm— — cos— ;— , sm— — sm— — , cos— > (7.8) 
k k k k k J \ ! u=o,....k ^ ' 



=0,...,fc-l 



of which the total number is n = fc^ — /c + 2. 



One may be overwhelmed by the number of the elements of 1^ because solving the (fc^ — fc + 2)- 
dimensional problems is required here. There is, however, a way to solve this part of problem. The 
calculation of H7.7|l is performed by the mathematical programming package NUOPT developed by 
Mathematical IncESnExl. NUOPT is capable of solving large optimization problems with an interior 
point method. The computation time is approximately one hour for n = 1562 (fc ~ 40) and approximately 
one day for n = 4831 (fc = 70) on a Sun Microsystems workstation with a CPU UhraSPARC-II 360MHz 
and 2048 Megabytes of RAM. It seems to take time proportional to n^, or n to the power of three. 

As fc becomes larger, the finer the lattice becomes, and the smaller the coarseness 8 becomes due to 
(5 — 0(kr^). Fig. 17.41 Right shows how each solution of the approximated problems 1)7.7(1 converges to 
the Holevo capacity of the qubit channel 

A4 : p(x, 2/, z) ^ p{0.6x + 0.021, 0.601?/, 0.5^ + 0.495), (7.9) 

which is a distinct one, as the author calculated the Holevo capacity with quite high precision, of which 
the detail is described in Section 1731 

7.2.3 Error analysis 

The following theorem is significant in that it guarantees the convergence of the approximated 
solution to the actual Holevo capacity, i.e. c^^ C{A) as fc ^ 00. 

Theorem 7.1 

The error, or the difference, between the Holevo capacity C(A) and the solution of the approximated 
solution of i7.7}) is upper-bounded by 0{—S^ log (5), which is equal to 0{k^^ logfc) when the lattice is £k- 
Particularly, when the output ellipsoid does not touch the surface of the Bloch sphere, it is bounded by 
0(-(52), which is 0{k-^) for Ik- 
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y=0.05/kA2 




Left: A lattice ^20 on the Bloch sphere Left: How the calculated values converges 

Figure 7.4: Left: A lattice covers dB{Tiout)- To solve the maximization problem (|7.7|l . {pi\i is specified 

beforehand, such as i?20- (Other lattices such as and £100 appear in Fig. 12.31 ^ 

„ r / the concavity w.r.t. (pi)i of the objective function, \ ., . „ -11, • • ,1 i • 

Because of ( 1 1 1 \ >, it is feasible to maximize the obiective 

\ and the convexity 01 the search space (pi)i ' 

function. One can get a value close to the actual Holevo capacity, when the coarseness on the lattice is 
small enough. 

Right: The actual convergence for an example channel. The horizontal axis is a log plot of k, and the 
vertical axis is a log plot representing the difference between the calculated value and the optimum value. 
A line y = 0.05/fc^ is drawn for reference. 



Proof 

Denote Ck as the solution of the approximated problem i)7. 7|) with the lattice £k, and Cao as the 
Holevo capacity C(A). Then the following conditions hold: 



Ck < Coo 

Ck = max{H (Ap 1 1 Aak) ; p € ik} 
Coo < max{H {Ap\\ Aak) ; p € dBiHin)} 

where at = '^iPiPi, the average of calculated engaging inputs of ^7. 7|) with Ik- Denote 

fk{p)^H{Ap\\Aak) for p^OBiHt.). 

The above equality and inequalities i)7. JCIj) . i)7.ii|) and \7.1'J^ lead to 

\coo- Ck\ < max fk{p) - max f kip). 



(7.10) 
(7.11) 
(7.12) 

(7.13) 

(7.14) 



JVote that £k C dB{Ti.±n) and the coarseness 6 is the upper bound of the distance from any points 
on dB{Hin) to the nearest point among £k- Therefore, estimating RHS of 1^7.14}) is reduced to: 



the infinitesimal analysis of Jk{') with the domain dB{TL±-a) around the maximum. 
We need to consider the following two cases: 



(7.15) 
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1. When the maximum is given at pq where ApQ e B{Hout)\dB{Hout), /(■) is as doubly continu- 
ously differentiable as _ff (• || Aak) on B{Ti.out)\dB{Hout)- Note that the compactness ofdB{Ti.in) 
bounds the coefEcient of the second derivative of f{-). These lead to \7.14\ . RHS) — 0{5^). 

2. Otherwise (when the maximum is given at 9S(7iout)J, one needs to consider the infinitesimal 
behavior of H {■ \ \ Aak) around dB{TLont)- This infinitesimal behavior is as same as the infinites- 
imal behavior of the function — zlogz substituted by z = around t = Q. This leads to \7.14\ . 
RHS) = 0{-6^ log (52) = 0{-S^ log (5). 

Combining these two results above, RHS of 1^7.1^) is upper-bounded by max{0{6^),0{~6^ \ogS)} = 
0{-S^ log 6). Namely, 

\c^ -Ck\^0{-5'' logS) (7.16) 
Since 6 = 0{k-^), \coo ~ Ck\ = 0{k-^ logfc) holds. 




□ 

7.3 Application — finding qubit channels with four engaging inputs 

The channel 

A4 : p{x, y, z) ^ p(0.6a; + 0.021, 0.601y, 0.5z + 0.495), (7.17) 

already appeared in (|7.9|l has a distinguished property: The number of its engaging inputs is four. As 
far as the author checked tens of times by the method presented in the previous section, the number 
of the engaging inputs of qubit channels are, mostly two, rarely three. Four does not appear without 
elaborate preparation. 

7.3.1 When a qubit channel has four engaging inputs 

When all the approximated positions of the engaging inputs of a given qubit channel A are located 
through the method described in the previous section, one can calculate its Holevo capacity C(A) much 
more precisely by the Newton-Raphson method, even though the significant precision seems limited to 5 
- 6 digits by the method in the previous section. We employed this method for the channel A4, and the 
capacity and the engaging inputs are shown in Table. ITtI The calculation is performed by mathematical 
software Maple 6, and we actually obtained the precision of C{AiC] as much as more than 50 digits. 

Incidentally, by Caratheodory's theorem, the theoretical upper bound of the engaging numbers of 
qubit channels is four jHayashi-imai-Matsumoto-Huskai-fcihimono U4| ^ Therefore, the discovery of this channel 
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Kip{x, y, z) = p(0.6a: + 0.021, 0.601?/, 0.5z + 0.495) capacity = 0.3214851589 


probability 


engaging input {x,y,z) 


e 


0.2322825705 ( 


0.2530759862, -0.0000000000, 0.9674464043) 


14.66° 0.00° 


0.2133220819 ( 


0.9783950999, 0.0000000000, 0.2067438718) 


78.07° 0.00° 


0.2771976738 {- 


-0.4734087533, 0.8646461389, -0.1681404376) 


99.68° 118.67° 


0.2771976738 {- 


0.4734087533, -0.8646461389, -0.1681404376) 


99.68° -118.67° 


average { 


0.0050428099, 0.0000000000, 0.1756076944) 






0, 6 denote the angular coordinates of the engaging inputs. 


probability 


engaging output {x, j/, z) 


5[A4(p)] 


0.2322825705 


( 0.1728455917, 0.0000000000, 0.9787232022) 


0.0300135405 


0.2133220819 


{ 0.6080370599, 0.0000000000, 0.5983719359) 


0.3786915585 


0.2771976738 


(-0.2630452520, 0.5196523295, 0.4109297812) 


0.5935800377 


0.2771976738 


(-0.2630452520, -0.5196523295, 0.4109297812) 


0.5935800377 


average 


( 0.0240256859, 0.0000000000, 0.5828038472) 


0.7383180644 


Table 7.1: Data for the qubit channel A4 with the engaging number four 




Figure 7.6: The quantum divergence H (A4P || A40') on the surface of the output ellipsoid of A4. p runs 
on the surface of the Bloch sphere and the a is the optimal average input. The four most reddish points 
correspond to the engaging output because they are the farthest points from K/^a w.r.t. the quantum 
divergence. The rectangular region is a projection from the ellipsoid surface. (The ellipsoid in the figure 
is drawn slightly transparent to see through to the opposite side.) 
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A4 settle the problem of the actual maximum engaging number for qubit channels. It is the following 
discovery after jKmg-Nathanson-Ruskai 02| -^yjjjcjj found qubit channels which require three inputs to achieve 
the Holevo capacity. Moreover, it is a kind of the final goal that jF^chs 97 | attempted to dispel a prejudice 
that the engaging inputs should be orthogonal states, i.e. two antipodal points on the Bloch sphere; 
the idea is that a general quantum channel is noisy, therefore the input should be kept to the most 
distinguishable states. On that prejudice, for any qubit channel, the engaging inputs should be two. 

One of the reasons that non-orthogonal inputs achieve greater capacity is that collective measurement 
on the output is allowed. This consideration may lead to the idea that the collective preparation of inputs, 
which are entangled input states, may achieve further capacity than the Holevo capacity. This leads to 
the additivity check of the Holevo capacity, described in Subsection l7.3.3l 

7.3.2 How w^as A4 found? 

How the channel A4 : p{x,y,z) 1-^ p(0. 6a: + 0.021, 0.601j/, 0.5z + 0.495), which has four engaging 
inputs, was found is explained here. 

- 1. The author perceived the tendency that the engaging outputs appear near the two endpoints 
of the major axis of the output ellipsoid. Moreover, when the surface of the output ellipsoid 
very closely approaches the surface of the Bloch sphere, then the third engaging output on the 
ellipsoid may appear near the Bloch sphere. 




# = 2 # = 3 



- 2. The output ellipsoid of an already-found three-engaging channel A3 : p{x, y, z) i— s- p(0.6x, 0.6y, 

0.5z -I- 0.5) is stretched along the y-direction very slightly to an adequate extent. The channel 
became A3' : p{x, y, z) ^ p(0.6x, O.GOly, 0.50-1-0.5). It is desired that one of the engaging out- 
puts splits into two (we denote this phenomenon bifurcation)] however, the excessive symmetry 
in terms of the location of the ellipsoid prevents bifurcation. 

- 3. To reduce the symmetry, moving ahead in the x-direction is attempted, and to keep the output 

ellipsoid inside the Bloch sphere, one must pull the ellipsoid to the center slightly along the 
z-direction. The channel became Ky^ ■ p{x, y, z) 1-^ p{0.6x + e, 0.601?/, 0.5z -I- 0.495). 

- 4. Various e were tried, and bifurcation occurred when e = 0.021. 




# = 3 # = 4 
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- 5. CPTP-ness was checked by the method introduced by i R"«kai-Szarek-Werner 02 |^ j,^^ method 

showed |e| < 0.05277.. is equivalent to Aye being CPTP jHayashi-Imai-Matsumoto-Ruskai-Shimono 04| ^ 

Another qubit channel with four engaging inputs: 

A : p{x, y, z) ^ p{0.8x + 0.022, 0.8015y, 0.750 + 0.245) (7.18) 
was found by the prescription above from a channel A : p{x, y, z) i— > p(0.8x, 0.8?/, 0.75z + 0.25). 

7.3.3 Additivity check for A4 (g) A4 

As mentioned in Subsection 17.3.11 a quantum channel A can get more communication efhciency if 
collective measurement by the receiver is allowed. This efficiency is called the Holevo capacity of this 
channel. One might like to consider whether "collective preparation" by the sender increases further 
efficiency, i.e. whether input particles with quantum entanglement increase the efficiency. This is for- 
malized as whether C(A'^") > nC{K) for a channel A and a positive integer n. We checked this for 
n — 2 with A4, which is an eccentric channel as already stated. 

There is a useful relation to check the additivity iHayashi-lmai-Matsumoto-Huskai-Shimono 04| 

C{K®^) > 2C(A) ^ max H (h®^uj \\ > 2C(A), (7.19) 

where cr^ is the optimal average output of A. It is useful in that it greatly reduces the dimension of the 
search space. Therefore, lu = is substituted by 



1^) = ^\u)(g, |f) +e^Vl-p|"^> ® \v^) (7.20) 

-here |.) = (,^ll,J» = (.^ll.J, l"^> = {'-Zllt'^) \v-) = {^^^t'^)- 
( are orthogonal vectors of \u) and \v), respectively.) 

The range of the parameters of H7.20|l are 

< p < 1, < 9u,0y, iy<2Tr,0< 0„ < |. (7.21) 

(|7.20|) is simply an ordinal Schmidt decomposition form of a bipartite pure state. We plot H (^An'^'^uj \ \ c 
with various u and w as a function of p and v. A typical plot is shown in Fig. 17.71 All the graphs of 
the numerical experiments show that i7 ^A4^^cj 1 1 cr^^'*^^ is deeply convex w.r.t. p, the entanglement 
parameter of the Schmidt decomposition. Therefore, the extra gain in the Holevo capacity of A4^^ does 
not seem expected. 

(As a result, such deep convexities may lead to whether this convexity universally holds for any 
quantum channel. If it holds, the additivity holds. However, this convexity does not always hold. The 
convexity in H {p \ \ a) as of p leads to the concavity of S{p); however, as Fig. l7.8l shows. S{A'^'^\tp){ip\) of 
A : p{x,y,z) p(0.75x, 0.75?/, 0.5z) with = .^|00) + -y/1 — p |11) is convex as of p. Therefore, the 
additivity problem is not yet solved. ) 

7.4 Conclusion and discussion 

7.4.1 The calculating method vifith universal applicability and convergence 

The calculations of the Holevo capacity were performed by the author for the first time in the world 
in the following sense: 
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If there is an extraordinary gain in the communication 
capacity of A4^^, there exists a pure state cu shared 
by two particles which satisfies H ^A^^'^uj \ \ > 
2C{A/i). Our numerical experiments, however, show 
that the more the two particles are entangled, the 
more the communication capacity falls, as the figure 
shows. 



Figure 7.7: Convexity of H (^A^"^^ uj \\ a'^^^®^^ 



One might like to conjecture that entanglement al- 
ways causes convexity as illustrated in Fig. 17.71 
Note that the convexity of H{A^'^Lo\\a'/^^^) causes 
the concavity of S{A'^'^uj) because H{A'^^uj\\a'^'^^) + 
5'(A®^aj) is linear as of co. The figure on the left, 
however, shows the opposite for A : p{x,y,z) 
p(0.752;, 0.75?/, 0.5^) with It/)) = y^|00) + yr^|ll), 
as of a function of p. 



p 

Figure 7.8: A counterexample of the concavity conjecture 



- 1. The calculation method was applicable to any qubit channel. 

- 2. The answer was truly close to the real answer; convergence to the real answer was guaranteed. 

The method involved a high-dimensional optimization and was actually performed for approximately a 
hundred channels. 

The convergence speed is, however, not so fast in that bounding the error within e > requires a 
lattice with the coarseness of 0(e^/^+°). If this lattice is £k, k — il(e^^^'^), thus n — 8(e^^^") where n 
is the number of the points of ik- Assuming the calculation cost in time is Q{n^), it ends up to become 
equal to 0(e^'^^''). (Here +0 and —0 mean arbitrary positive and negative numbers, respectively, with 
the intention that they are non-zero and are close to zero.) 

7.4.2 An attempt to check the additivity in Holevo capacity 

To investigate a mysterious aspect of quantum physics, the additivity of the Holevo capacity was 
checked for a channel with four engaging inputs, which was discovered by the proposed calculation 
method. Since the communication efhciency through a quantum channel increases by collective mea- 
surement on the receiver side, one might expect it to also increase by the collective preparation, i.e. 
entanglement among inputs may increase the efficiency. The discovered channel was such a candidate, 
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as it increases efficiency by the collective measurement by the receiver, but the analysis presented in this 
chapter showed the collective preparation of the sender cannot seem to increase the efficiency. 
In order to refine this attempts, we may need to investigate: 

- 1. how the input and the output of a quantum channel are reversed, and 

- 2. how communication efficiency relates to the coUectiveness of both the input and the output of 

a quantum channel. 

7.4.3 Beyond qubit channels 

One may like to calculate the Holevo capacity of a quantum channel which is not a qubit channel. 
Because of dimensional explosion, it is impossible to calculate the Holevo capacity even for three-level 
systems. The reason is as follows: the corresponding space of the Bloch sphere for c?-level system forms 
a {(P — l)-diniensional manifold which is the convex hull of a {2d — 2)-dimensional manifold. What 
will happen if we cover the {2d — 2)-dimensional manifold with a mesh of 40 x • ■ ■ x 40 ? A roughly 

approximated number of vertices for each d is as shown in Table 17.21 Thus it is quite hard to calculate 



d 


2 


3 


4 


8 


16 


2d -2 


2 


4 


6 


14 


30 


4Q2d-2 


1.6 X 10^ 


2.6 X 10*^ 


4.1 X 10^ 


2.7 X 10^2 


1.2 X 10^** 



Table 7.2: Vertices number of the lattices for each quantum channel of d-level system. 

the Holevo capacity beyond qubit channels. A quite elaborate methodology may be required to realize 
this calculating, for example, an efficient configuration of the vertices of the mesh. 

Endnotes 

1 Classical information capacity: it is classical, the antonym of quantum. The Holevo capacity 
measures the efficiency of transmitting the sequences of the alphabets of some message from the sender 
to the receiver through a given quantum channel, while the quantum capacity may be defined as the 
efficiency of transmitting the input quantum states themselves with high fidelity through the channel. 




Transmitting Classical information Transmitting quantum information 
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Part IV 

Epilogue 



Chapter 8 



Conclusion 



8.1 Summary of this dissertation 

In this dissertation, the following concepts in quantum information science have been dealt with: 

- 1. the entanglement cost and 

- 2. the Holevo capacity. 

These two concepts are important quantification in the emerging science of quantum information. To be 
slightly more precise, they are both quite important candidates to measure 

- 1. the anomalous correlation between multiple figures peculiar to quantum physics, and 

- 2. the ultimate communication efhciency of given channels in this physical world, 

respectively. Note that the anomalous correlation, which is called quantum entanglement, is an 
indispensable resource in quantum computation and quantum communication. Also note that the two 
concepts are related theoretically through Stinespring's dilationlMi i""'°t°-Shimono-Wfinter 04 | _ 

From an objective overview, the desirable goal in investigating entanglement measures and capacity 
of quantum channels might be as follows: 

- identifying the most approvable entanglement measure if it exists, 

- discovering how to calculate this measure, 

- determining whether measures based on various concepts are equal to each other, 

- determining the additivity of the measures (perhaps e{p<S)(T) 7^ e(p) + e(cr) would occur for some 
measure e contrary to expectation), and 

- determining whether the Holevo capacity is the most approvable communication capacity of 
quantum channels; this issue is equivalent to whether the additivity of the Holevo capacity holds. 

- discovering how to calculate the communication efficiency of quantum channels. 

With the purpose of settling the problems written above, this dissertation has attempted to solve 
the following problems: 

- Ec (the entanglement cost) and Ep (the entanglement of formation) of antisymmetric states, 

- the difference between Ec and i?£)(the entanglement distillation), 

- the superadditivity of Ep, 

- the algorithm to calculate the Holevo capacity, and 

- the additivity of the Holevo capacity. 
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8.2 Application in daily life 

— A long, long time from now or in a galaxy far, far away — 

How will the research of this dissertation contribute to the real world? One of the amusing applica- 
tions is quantum teleportation |sHS^^^ 



Quantum Entanglement 




A Cricket Teleports ! 

Figure 8.1: A sketch of quantum teleportation. 

A quantum state (a cricket in the sketch) in Alice's hand is teleported to Bob with quantum entanglement and 
classical communication. When the age of quantum teleportation comes, each of the amount of the quantum 
entanglement and the classical communication in quantum teleportation is important. 



Quantum teleportation is a method to teleport the quantum state in one place to a remote place 
with quantum entanglement shared by the two sites beforehand. You can consider Alice desires to send a 
quantum state to Bob in a remote place. They need to have pairs of entangled particles divided between 
Alice's site and Bob's site. The quantum state to be teleported may be anything, but it is desirable to 
be a small object like DNA molecules or mineral samples on asteroids. Alice may teleport herself with 
a huge amount of entanglement and with a transmitter of huge information. A brief explanation of the 
method to teleport the object is as follows (See Fig. I8.1|l : 

Alice makes the object and the particles in her hand interact with each other. After getting 
signals from the interaction, she sends the signals to Bob. According to the signals, he 
manipulates the particles in his hand which turn into the object Alice once had. 



How are these operation performed? Here is a theoretical description how to teleport the object 

)o ~ Si=o '^i N)o' 



of d-level system of which the basis is a {\i)}i=o.i,...,d- Assume the state of is \il')o — SiLo^ '^i N)o' ^^'^ 



Wi}'iZQ is not need to be known. In the following, ui = e , and + operation is done in modulo d. 

- 1. The shared entangled particles AB are prepared in X^it^o^ N)a ® N)b- 

- 2. The interaction is done by 'measuring' the particles in her hand to physically identify with one of 



\''Px,y) = —/= \i)o |i + y)f^ > ; {\ipx,y)} is an orthogonal basis of the space of OA. 

y=0,...,d-l 



3. Bob manipulates his particles B with a unitary operator Uxy — (uij) where = Si^+y lo^^ and S 
is the Kronecker delta (Skk = 1 and S^t = when k ^l). 
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(For d = 2, the entangled state shared beforehand is the Bell state, the measurement is the Bell mea- 
surement, and the unitary operator is the Pauli operator.) 

When Alice measures that the state of OA is \ipx.y)Q/:^, the state of the particles in Bob's hand turns 
to be \(j>x,y)B = Efc =0 " '^i I* + ?/)b' because the state of the whole OAB before the measurement is 

^ -j= |V-'t,j/)oa ® I't'x.y)^- He retrieves the state of X]f=o^ '^i N)b from \(f)x.y)Q by the unitary operation 
with the information x and y. When the physical attributes, such as atomic configuration or the state 
of electrons, of \i)Q and |z)g correspond, one observes the object in Alice's site teleports to Bob's site. 

The technical difhcultics are as follows: how to cope with the decay of the entanglement (except 
photons, the quantum state is quite easily interacts with its outer environment within a millionth second, 
which causes quite hard to preserve quantum entanglement), and how to transmit the huge information 
(even a Imm^ (cubic millimeter) of water requires d > 10^°^°, thus Alice needs to tell x and y to Bob 
with huge digits. It takes much more than a trillion years by the present technology). Finding the 
technology to preserve quantum entanglement for a long time and to load information on each light 
particle of ordinary brightness will solve these difficulties. 

When the age of quantum teleportation comes, measuring quantum entanglement and also the 
ultimate capacity of information- transmitting channels will become important. They can be compared 
to the amount of petroleum and the road capacity for an automobile. The research of this dissertation 
will hopefully contribute to the furthering understanding these points. 
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= Tr {p log p - p log a) 

(l:^Kullback-Leibler£ggt) 
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7T(A®2|^)(^| II (Ap)®2) 
\^) = ^/V\u)® \v) 
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